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overall  performance  in  systems  that  respond  at  vibration  frequencies.  Stick 
location  had  no  significant  effect  on  either  tracking  or  biodynamic  performance 
measures.  The  vibration -correlated  component  of  tracking  error  was  relatively 
small,  vibration  effects  were  accounted  for  largely  by  changes  /.n  pilot  model 
parameters  related  to  visual  resolution,  time  delay,  and  motor-, -elated  remnant. 

A model  based  guide  for  the  design  of  control  sticks  in  a vibration  environment 
is  described.  This  model  is  based  on  the  state-variable  model  for  pilot  /vehicle 
systems.  Effects  of  vibration  are  represented  as  additional  model  element, s and 
by  changes  in  pilot-related  parameters  of  the  tracking  model. 
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The  experiments  W(jre  conducted  at  the  vibration  branch  of 
the  Aerospace  Jlcdical  research  I.aboratorv.  They  v>»ere  made 
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f^r.  Charles  Clausor  of  t.he  AMRL  kindly  took  the  anthropometric 
measurements  of  the  suy>;ects. 

Lt,  Philip  D.  Houck  was  the  technical  monitor  for  this 
contract  and  project  engineer  for  the  AMRL  experimental 
program.  Lt.  Houck  is  the  author  of  Section  2 of  this 
report. 
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1.  INTRODUCTION  AND  SUMMARY 


1.1  Objective  of  this  Report 

In  recent  years  aircraft  have  been  constructed  such  that 
they  can  perform  in  environments  that  severely  tax  or  exceed  the 
limits  of  human  tolerance.  Fevere  vibration  is  one  such  environ- 
ment. 


Whole-body  vibration  acts  in  a number  of  ways  to  cause 
performance  degradation  of  manned  vehicle  systems.  Vibration 
is  transmitted  to  the  seat  of  the  pilot  by  the  aircraft  and 
propagates  through  the  pilot's  torso  and  arm  to  the  control 
stick  to  produce  control  inputs  that  are  linearly  correlated 
with  the  vibration  input.  Vibration  also  increases  the 
stochastic  portion  of  the  pilot's  control  activity  (i.e., 
"remnant")  due  partly  to  visual  interference  effects  (e.g,, 
blurring  due  to  relative  eye-display  motion)  and,  to  a greater 
extent,  to  motor  interference  effects  such  as  noise  injected 
into  proprioceptive  feedback  paths. 

The  effects  of  vibration  on  control  performance  are  to 
some  extent  influenced  by  the  design  of  the  control  stick. 
Accordingly,  the  Aerospace  Medical  Research  Laboratory  has 
recently  completed  a series  of  man-in-the-loop  experiments  to 
determine  mathematically  the  relation  between  control-stick 
design  parameters  and  tracking  performance  in  a vibration 
environment.  The  objective  of  this  report  is  to  summarize  the 
results  of  this  experimental  study  and  to  propose  guidelines 
for  the  design  of  control  sticks  in  vibration  environments. 
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1,2  Summary  of  Experimental  Results 

The  primary  variables  of  the  experimental  program  were  (a) 
stick  desian  parameters,  (h)  stick  location,  and  (c)  presence  or 
absence  of  whole-body  vibration.  Six  stick  configurations  were 
explored;  three  "sprino  sticks"  havina  relatively  low  spring 
gradients,  and  three  "stiff  sticks"  havino  suljstantially  larger 
spring  qradients.  Side  and  center  stick  locations  were  explored. 

single  z-axis  sum-of-sines  vibration  environment  was  used  in 
this  program,  and  the  biodynamic  environment  (pilot/seat  interface, 
seat/platform  interface,  etc.)  v;as  kept  as  invariant  as  possible. 
Some  of  the  important  experim.ental  results  are  summarized  below. 

I' f feats  of  Stick  Parameters , Stick  design  had  relati- 
vely little  effect  on  tracking  error  under  vibration 
conditions  for  the  aircraft  dynamics  that  were  used, 
but  considerably  greater  control  activit^'  was  observed 
for  the  stiff  sticks  than  for  the  spring  sticks.  This 
control  activity  could  cause  problems  if  the  system  or 
sul:)system  responds  to  those  frequencies. 

Vibration  Feedthrouah . Stick  dcsian  naresmeters  influenced 
primarily  tlic  vil)ration-correlated  (i.e.,  "feedthrough") 
component  of  tracking  error.  These  effects  can  l)e 
represented  Ijv  an  impedance  model  v;hich  includes  the 
stick  Impedance  plus  tv;o  imnedanco  functions,  inde- 
pendent of  stick  design,  that  account  for  biodvnamic 
response  behavior. 

Vibration  Interference  Effects.  Feedthrough  accounted  for 
a small  fraction  (less  than  10?.)  of  the  tracking  error 
variance.  Tn  terms  of  the  pilot/vehicle  model  emplovcd  in 
tlie  analysis,  the  primary  effects  appear  be  an  increase 
in  motor- related  remnant  and  an  increase  in  time  delay. 

Visual  effects  appear  to  be  secondary. 


Model  Results . Good  agreement  between  model  outputs 
and  experimental  measures  were  obtained  for  the  six 
control  stick  configurations  with  a fixed  set  of 
model  parameter  values. 

Stick  Location,  Stick  location  had  no  significant  effect 
on  rms  tracking  score,  nor  were  significant  effects 
observed  on  rms  head,  shoulder,  or  elbow  accelerations. 

Individual  Differences . Individual  differences  with 
regard  to  biodynamic  response  measures,  relative  to  mean 
performance,  were  generally  no  greater  than  subject-to- 
subject  differences  in  tracking  performance.  Furthermore, 
not  all  such  differences  were  statistically  significant. 
Thus,  it  is  reasonable  to  analyze  performance  of  the 
"average  raan", 

1.3  Nature  of  the  Deslcin  Guide 

If  we  were  to  follow  traditional  format,  the  design  guide 
would  consist  of  a set  of  tables  and  charts  to  allow  rapid, 
straightforward  optimization  of  control  stick  parameters  for  a 
variety  of  aircraft  and  for  a variety  of  vibration  and  biodynamic 
environments.  Perhaps  the  greatest  obstacle  to  the  development 
of  a guide  of  this  sort  is  the  inability  to  extend  existing 
results  to  other  vibration  and  biodynamic  environments. 
Comprehensive  measures  of  tracking  performance  have  not  been 
obtained  in  a wide  variety  of  such  environments,  and  reliable 
theoretical  models  for  biodynamic  response  behavior  and  for 
vibration-related  interference  effects  are  lacking  to  allow 
extrapolation  cf  existing  results. 


Another  potential  problem  in  posed  by  the  complexity  of 
the  interaction  between  control  task  parameters,  stick  design, 
vibration  parameters,  and  biodynamic  response  mechanisms.  All 
these  factors  interact  in  ways  that  are  not  entirely  linear  to 
influence  tracking  performance.  Thus,  it  is  not  clear  at  this 
point  that  a handbook-type  guide  could  be  provided  for  general 
application  without  requiring  an  inordinate  cataloguing  of  task 
situations. 

We  therefore  propose  in  this  document  a design  guide  based 
on  a rather  extensive  model  v’ith  dem.onstrable  predictive  capabi- 
lities. With  such  a model,  the  designer  can  not  only  optimize 
the  design  of  the  control  stick,  but  he  can  also  explore  the 
effects  of  other  important  factors  on  pilot  behavior  and 
overall  system  performance. 

The  model-based  guide  proposed  in  this  report  is  built 
around  the  state-variable , or  “optimal- control"  model  for  pilot,' 
vehicle  systems.  Vibration  effects  are  accounted  for  by  addi- 
tional model  elements  and  by  changes  in  pilot-related  parameters 
of  the  tracking  model.  The  procedure  f'*)r  using  the  proposed 
model  to  predict  t,'.ie  effects  of  stick  parameters  on  tracking 
performance  is  diagrammed  in  Figure  ''  ind  io  suinmariitd  below: 

1.  System  dynamics  are  represented  in  state-variable 
format.  Control-stick  and  display  dynamics,  as  well 
as  frequency-shaping  of  the  tracking  input,  are 
included  in  this  formulation. 

2.  Pilot-related  model  parameters  not  affected  by 
vibration  are  assioned  values  from  v/ell-documented 
studies  of  pilot/vehicle  performance  in  non- 
vibration environments. 
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riaure  1,  Procedure  for  Predicting  Tracking  Performance 


f:  .1  ■ 
I' ' '. 


3.  Model  parameters  relating  to  pilot  time  delay 

and  motor-related  sources  of  pilot  remnant  are  assianed 
values  to  reflect  the  effects  of  the  vibration, 

4.  Describing  functions  relating  head  motion  to 
platform  acceleration  are  used  to  determine 
relative  motion  between  the  eye-point-of-regard 
and  the  display.  This  motion  variable  determines 
the  value  of  a model  parameter  related  to  visual 
resolution  limitations. 

5.  Having  quantified  all  parameters  of  the  pilot/vehicle 
model,  predictions  of  tracking  V'^erformance  (including 
remnant)  are  nov;  obtained. 


i.< 


6.  The  impedance  model  for  stick  feedthrough  is  analyzed 
to  predict  vibration-correlated  stick  motion  for  the 
stick  configuration  of  interest,  and  the  variance  of 
the  tracking  error  due  to  feedthrough  is  computed. 

7.  Error  variance  related  to  feedthrough  is  adde'd  to 
the  variance  score  with  tlie  pi lot/vchic le  model  to 
yield  total  error  variance. 

8.  Since  a change  in  stick  parameters  will  generally 
affect  both  tracking  performance  as  well  as  feedthrough, 
the  entire  procedure  is  repeated  to  explore  the  effects 
of  changing  one  or  more  parameters  of  the  control  stick. 


■i 
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Once  reliable  theoretical  models  for  blodynamlc  response 
mechanisms  ana  vibration-related  Interference  effects  become 
available/  It  will  be  possible  to  use  this  model  to  explore 
the  Interactions  between  optimal  stick  design  and  factors 
such  as  vehicle  dynamics/  vibration  parameters/  and  blodynamlc 
environment.  If  some  of  these  Interactions  are  found  to  be 
unimportant/  or  If  they  can  be  accounted  for  In  a simple  manner/ 
It  may  then  be  possible  to  define  a set  of  guidelines  In  handbook 
format. 

1.4  Limitations  of  the  Design  Guide 

Because  we  lack  reliable  models  for  extrapolating  existing 
results  to  other  vibration  and  blodynamlc  environments/  the 
design  guide  described  In  this  report  Is  intended  to  apply 
strictly  to  the  biodynamic/vibration  environment  employed  in 
the  recent  experimental  study.  Because  of  the  demonstrated 
capability  of  the  pllot/vehlcle  model  to  predict  tracking 
performance  in  a variety  of  control  situations/  however/  we 
expect  the  guide  to  be  useful  In  exploring  tracking  tasks  and 
stick  designs  beyond  those  studied  In  the  laboratory.  If  the 
guide  is  to  be  extended  to  other  blodynamlc  and/or  vibration 
environments/  "calibration^’  experiments  will  have  to  be  per- 
formed to  determine  the  effects  of  vibration  on  pilot  parameters 
and  to  quantify  relevant-  b’odynamic  response  mechanisms. 

The  reader  should  be  aware  that  the  control-Mtick  design 
influences  primarily  vibration  feedthrough/  which  typically 
accounts  for  a small  portion  of  the  increment  in  tracking 
error  due  to  vibration.  Stick  design  is  an  important  consi- 
deration nonetheless/  for  an  inappropriate  choice  of  parameters 
can  result  in  much  greater  feedthrough  than  is  necessary  and/or 
require  control  forces  or  displacements  that  are  outside  the 
rang^  of  efficient  tracking  performance. 
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1.5  Suggested  Areas  for  Further  Research 

The  following  areas  of  research  are  suggested  for  improving 
our  understanding  of  the  relationship  between  tracking  performance 
and  vibration  and  thereby  allowing  developnent  of  a more  general 
design  guide. 

Biomeohaniaal  Models,  Considerable  effort  is  reouired 
to  develop  reliable  theoretical  models  fo]'  biomechanical 
response  mechanisms  so  that  the  effects  of  vibration  on 
head,  body,  and  stick  motions  can  be  predicted  for  a 
variety  of  biodynamic  and  vibration  environments. 

Tracking  Interference  Effects,  Theoretical  models  are 
desired  for  predicting  the  effects  of  vibration  on  changes 
in  pilot-related  model  parameters  ~ particularly  motor- 
related  remnant  parameters  and  time  delay.  Further 
experimentation  with  different  vibration  spectra  and 
different  biodynamic  environments  would  be  needed  to 
devcl'  n such  models. 

Further  Development  of  Pilot/Vehicle  Model,  The  existing 
state-variable  pilot/vehicle  model  is  in  a high  state  of 
development  in  that  performance  caa  be  predicted  for  a 
variety  of  tracking  situations  with  a consistent  set  of 
model  parameters.  Nevertheless,  further  theoretical 
study  is  suggested  to  refine  aspects  of  the  model 
related  to  motor  limitations;  specifically,  to  isolate 
one  or  more  parconeters  that  account  entirely  for  (and 
only  for)  degradation  of  motor-related  sensory  information. 
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I Biodynamio  Response  Linearity , Experimental  studies  are 

h 

I needed  to  determine  the  extent  to  which  biodynamic 

I response  mechanisms  are  linear.  The  existence  of  such 

i linearity  would  greatly  facilitate  extension  of  the 

I proposed  design  guide  to  additional  vibration  environ-* 

i 

$ ments.  Experimentation  with  various  vibration  spectra 

j and  amplitudes  are  required  to  explore  linearity. 

i 
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1.6  Organization  of  the  Report 

The  main  text  of  this  report  is  devoted  primarily  to 
definition  and  illustration  of  the  design  guide;  detailed 

i 

exposition  of  experimental  results  is  reserved  for  the 
Appendices. 

The  biodynamic/vibration  environmenc  to  which  the  design 
guide  applies  is  defined  in  Section  2,  and  other  aspects  of  the 
experimental  program  are  described.  The  model  structure 
underlying  the  design  guide  is  described  as  generally  as 
possible  in  Section  3.  Section  4 presents  the  design  guide 
in  more  quantitative  terms,  and  three  numerical  examples  arc 
presented.  The  guide  is  validated  against  data  obtained  in 
the  experimental  program,  and  design  curves  for  stick  gain  and 
stick  damping  are  presented  for  a specific  tracking  situation. 

Analysis  procedures  used  for  data  analysis  are  presented 
in  Appendix  A.  Detailed  discussions  of  experimental  results 
are  given  in  Appendices  B (biodynamics  response  measures)  and 
C (tracking  performance) . Some  details  of  the  model  analysis 
are  given  in  Appendix  D,  and  we  review  methods  for  differentia- 
ting among  potential  vibration-r-  lated  interference  mechanisms, 
/'nthropometric  measurements  of  the  test  subjects  are  tabulated 
in  /'ppendix  F. 
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2.  EXPERIMENTAL  DESIGN 


2.1  Purpose  of  These  Experiments 

The  pilot  controlling  an  aircraft  in  vibration  environment 
shov^'R  degradation  of  performance  due  to  noisy  proprioceptive  feed- 
bad:  and  direct  vibration  feedthrouqh  from  the  shoulder  to  the 
control  stick.  The  primary  interfacing  of  the  pilot  to  the 
aircraft  is  throuoh  his  motor  control  system  and  the  control 
stick.  The  stiffness  of  the  stick  v/ill  determine  hov;  much  of 
the  viliration  signal  is  allowed  to  enter  the  control  loop  and 
also  gives  the  pilot  proprioceptive  feedback  accordino  to  the 
"feel"  of  the  stick.  The  purpose  of  these  experiments  v;as  to 
find  how  pilot  aircraft  performance  is  effected  by  control  stick 
parameters  such  as  location#  sprino  oradient#  damping#  and  elec- 
trical gain.  The  objective  is  to  provide  the  aircraft  design 
community  with  control  stick  design  criteria  that  will  allow 
pilots  in  particular  aircrafts  to  perform  optimally  in  a 
vibration  environment. 

2.2  Experiment  Design 

The  experiment  that  was  conducted  consisted  of  subjects 
flvino  a sinole  axis  compensatory  trackinc  task  with  a simple 
intoprator  as.  the  aircraft  dvnamics.  The  parameters  of  the 
experiment  wore  whole-body  vibration  or  ahsense  of  vibration 
and  six  control  sticks  varying  from  2 Ib/in  to  600  Ib/in  in 
spring  gradient.  PacJi  of  the  control  sticks  was  arranged 
in  both  the  side  and  center  location.  The  control  stick  gain 
for  each  of  the  six  sticks  v;as  hold  constant  and  is  listed  in 
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Table  4 (p.  54  ) . "’he  three  spring  sticks  had  a constant 
electrical  displacement  aain  of  2.5  volts/in.  Both  the  Ij.ght 
and  medium  stiff  sticks  had  the  same  electrical  force  gain 
of  ,9  volts/lb;  v/hereas,  the  heavy  stiff  stick  had  a gain  of 
only  .5  volt/lb. 

The  in  subjects  partici'pating  in  the  experiment  v;ore 
members  of  the  ha2ardous  duty  panel.  Seven  of  the  subjects 
participated  in  all  of  the  experimental  conditions  '.nd  these 
su})jects  wore  used  as  the  data  base.  The  body  structure  of 
the  subjects  varied  from  large  to  small  and  cover(*d  most  of 
the  range  of  present  aircrew*-  . Since  body  size  could  influence 
transmission  of  the  vibration  oignal  to  the  control  (e.g,,  a 
short  man  requires  areater  extension  of  his  arm  to  reach  the 
controls),  the  control  confiauration  v.fas  fixed  to  allow  inves- 
tigation of  the  effects  of  body  sire,  r table  of  anthropometric 
measurements  obtained  from  the  test  subjects  : s given  in  Appendix  B 

The  follov;ina  fiaures  and  diagrams  describe  the  experimental 
set  up.  The  trackincr  task  is  diagramed  in  r'iaure  2 and  the 
control  and  suJjject  orientation  on  the  FIXMODB  is  demonstrated 
in  Figures  2,  3-5.  The  bite  bar  that  v/as  held  between  the 
subjects'  teeth  was  used  to  detennine  rotational  and  transla- 
tional motion  of  the  head  sc  that  a corielation  could  be  made 
between  visual  dogradati^^n  and  head  motion  of  the  pilot,  /“n 
accelerometer  v;as  also  mounted  on  the  subject's  .ider  and 
elbov;  to  measure  the  vibration  transmitted  through  the  body. 

/'ll  data  ’/as  recorded  without  pre-processing  on  a 14-channel  FM 
recorder.  This  system  and  the  tracking  simulation  is  shown  in 
Figure  6. 
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The  vibration  environment  was  supplied  by  a large  amplitude 
mu)ti-degree  of  freedom  hydraulic  vibration  table  named  the  SIXM^de, 
It  has  six  degrees  of  motion  including  x,  y,  z,  roll,  pitch  and 
yaw.  This  table  has  a payload  capacity  of  1,000  lbs.  and  is 
capable  of  sinusoidal,  sum  of  sines  and  random  vibration.  The 
experiment  used  a sum  of  five  sines  vibration  with  equal 
acceleration  at  2 , 3,3,  5,  7 and  10  hertz.  The  RMS  level  of 
table  acceleration  was  approximately  0.3  g and  had  a crest  factor 
of  3 giving  the  subject  a maximum  peak  acceleration  of  about  1 g 
depending  on  hov/  the  sines  added  together.  The  acceleration 
input  is  very  clean  (for  a hydraulic  shaker)  because  of  the 
elastomeric  couplers,  which  filter  out  extraneous  hydraulic 
inputs . 

2.3  Vibration  Fee'^through  Example 

Vibration  feedthrough  and  its  possible  effects  on  control 
systems  is  demonstrated  in  the  time  traces  of  typical  error  and 
control  signals  shown  in  Figure  7.  Static  and  vibration  condi- 
tions arc  displayed  for  spring  sticks  and  stiff  sticks.  The 
vibration  feedthrough  is  nuite  evident  in  the  stiff  stick 
controller.  The  error  signal  contains  a much  smaller  feedthrough 
component  because  of  the  filtering  of  the  simulated  vehicle 
dynamics  (an  integrator);  nevertheless,  some  of  the  vibration 
peaks  can  hie  seen  in  the  error  signal.  With  the  spring  stick, 
the  dynamics  of  th*^  stick  itself  filters  much  of  the  vibration, 
giving  by  far  a -smoother  signal.  These  considerations  of 
vibration  feedthrough  are  very  important  and  sliculd  be  consi- 
dered in  tlie  deveJopment  of  all  new  systems. 
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Figure  7.  Time  Traces  of  Error  and  Control  Signals 


16 


2.4  Predi cti vfiness  of  the  Model 


In  addition  to  baseline  experiments,  another  was  performed 
to  test  the  predictiveness  of  the  model  that  wat  developed  in 
this  study.  The  results  of  the  experiment  and  prediction  are 
displayed  in  Table  1.  For  this  case,  the  "liaht  stiff"  stick 
confiauraiion  defined  in  Table  4 was  used,  with  the  electrical 
qain  reduced  from  0.9  to  0.45  volt/lb. 


Table  1 


Compariaon  of  Predicted  and  Measured  Performance 
"Liqht  Stiff*  Control  Stick 


m 


m 
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3.  MODEL  STRUCTURE 


3.1  Overview 

Because  tracking  performance  depends  on  the  specifics  of  the 
tracking  task,  the  properties  of  the  control  stick,  and  the  nature 
of  the  vibration  and  biodynamic  environments  in  a complex  manner, 
some  form  of  predictive  model  is  mandatory  if  we  are  to  analyti- 
cally explore  the  effects  of  these  parameters  on  performance.  We 
have  therefore  derived  a model-based  guide  for  the  design  of  control 
sticks  in  vibration  environments.  The  models  underlying  the  design 
guide  are  the  subject  of  this  section;  the  design  guide  itself  is 
described  in  Section  4. 

The  key  model  element  of  the  design  guide  is  a pilot/vehicle 
model  that  relates  system  performance  and  pilot  tracking  behavior 
to  elements  of  the  control  sy-"  tern  such  as  plant  dynamics,  control- 
stick  properties,  and  external  tracking  disturbances.  The  effects 
of  vibration  are  represented  partly  as  additional  model  components 
and  partly  as  changes  in  values  assigned  to  certain  parameters  of 
the  basic  pilot/vehicle  model. 

A skeletal  outline  of  the  model  structure  is  shown  in 
Figure  8.  For  simplicity  of  exposition,  we  consider  a control 
situation  of  the  type  explored  in  the  current  study;  namely,  a 
single-input,  single-output  control  system  subject  to  vibration 
disturbances  in  a tingle  axis.  This  model  can  be  readily  extended 
to  include  multi-input,  multi-output  systems. 

The  pilot  is  assumed  to  observe  a compensatory  display  of 
tracking  error  and  to  manipulate  a single  control  stick.  Because 
the  pilot  will  generally  extract  velocity  as  well  as  displacement 
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information  from  the  error  indicator,  tracking  error  is  shown  in 
Fig\ire  B as  a vector  quantity.  The  system  is  assumed  to  be 
disturbed  by  one  or  more  zero-mean  Gaussian  random  processes  "i" 
which  we  designate  as  "tracking  inputs"  to  differentiate  from 
vibration  inputs.  Tn  laboratory  situations,  tracking  inputs 
are  usually  added  in  parallel  with  the  pilot's  control  (to  simulate 
wind  gusts  acting  on  the  vehicle,  for  example)  or  in  parallel  with 
the  vehicle  output  (i.e.,  as  a command  input). 


We  consider  the  control  input  "6"  to  be  the  sum  of  two 
components:  5^,  the  control  input  that  is  intended  to  minimize 

tracking  error,  and  6^,  an  input  correlated  with  the  platform  vibra- 
tion input  Up  that  results  from  biomechanical  vibration  feedthrough. 
The  signal  6^  includes  both  the  response  to  the  tracking  error  as 
well  as  wide-band  stochastic  behavior  (i.e.,  "pilot  remnant"'. 

The  control  signal  "u"  that  is  actually  generated  by  the 
pilot  is  converted  to  the  electrical  input  *'(5^"  in  a manner  deter- 
mined by  the  mechanical  and  electrical  properties  of  the  control 
device.  If  the  electrical  output  of  the  device  is  proportional 
to  stick  displacement,  the  follov'?ing  relationship  obtains: 

6 (r)  = • U(s)/ZS(s)  where  is  the  constant  of  proportionality 

and  ZS(s)  is  the  forre/displacement  characteristic  of  the  control 
device  in  complex  notation.  Although  not  shown  explicitly  in 
Figure  3-1,  vibration  feedthrough  is  also  influenced  by  the 
properties  of  the  stick,  as  discussed  shortly. 

Pilot  randomness  is  represented  by  a set  of  observation 
and  motor  noise  processes.  The  pilot's  perceptions  of  error 
displacement  and  error  rate  are  assumed  to  be  perturbed  by 
white  noise  processes  v^,  and  an  additional  motor  noise  process 
v^  is  assumed  to  perturb  control  activity.  Other  pilot-related 
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limitations  include  an  effective  processing  delay  ("time  delay" 
jn  Figure  8 associated  for  convenience  with  sensory  inputs 
aid  a first-order  low-pass  filter  ("motor  lag")  applied  to  the 
pilot's  control  response.  Within  the  constraints  imposed  by  the 
pilot's  perception  of  the  task  and  by  the  limitations  mentioned 
above,  the  pilot  will  adapt  his  control  strategy  to  achieve  best 
{.erformance. 

Detailed  discussion  of  the  major  components  of  the  model 
follows.  In  Section  3.2  below,  we  review  methods  for  predicting 
pilot  behavior  and  system  performance,  and  we  describe  the 
rationale  for  quantifying  pilot  limitations  in  non-vibration 
situations.  The  status  of  the  biodynamic  model  elements  is 
discus jQd  in  Section  3.3,  and  we  discuss  modifications  to  the 
basic  pilot/vehicle  model  that  are  imposed  I'y  vibration  environments. 

3.2  Pllot/Vehlcle  Model 

Predictions  of  pilot  behavior  and  overall  system  performance 
are  obtained  with  the  aid  of  the  state-variable  (or  ''opt.imal‘’contrul" ) 
model  for  human  estimation  and  control.  The  model  is  baaed  on  the 
assumption  that  the  well-motivated,  well-trained  human  operator 
behaves  in  a near  optimal  manner  subject  to  his  inherent  constraints 
and  limitations.  Thus,  once  we  have  adequately  described  the  control 
system,  performance  requirements,  and  pilot  limitations,  we  can 
predict  pilot  behavior  and  system  performance. 

Detailed  mathematical  descriptions  of  the  model,  along  with 
supporting  experimental  data,  have  been  presented  in  the  litera- 
ture [1,  2].  Key  features  are  reviewed  below. 
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A linearized  description  of  system  dynamics  is  given  by  the 
following  state  equation: 

x(t)  = A x(t)  + B u(t)  E w(t) 


where  _x(t)  is  the  vector  that  describes  the  "state"  of  the  vehicle, 
u(t)  the  pilot's  (scalar)  control  input,  and  w(t)  a vector  of 
wide-band  ("white")  driving  noise  processes.  If  the  external 
disturbances  are  rational  noise  spectra  of  first  order  or  higher, 
the  resulting  "input  states"  are  included  in  the  state  vector  x(t) , 
as  are  states  that  may  be  needed  to  represent  dynamics  of  the 
control  device.  In  addition,  tracking  disturbances  introduced 
by  linear  mechanical  feedthrough  of  the  platform  vibration  can 
be  represented  in  this  manner  (provided,  of  course,  that  vibration 
feedthrough  can  bo  described  analytically) . 

The  display  variables  used  by  the  pilot  are  represented  as 
linear  combinations  of  state  variables  as  follows: 


y^(t)  = C x(t) 

wher>..  the  display  vector  ^{t)  typically  includes  both  displacement 
and  rate  variables  for  each  physical  display  indicator.  For  a 
single-variable  tracking  task  in  which  system  error  e(t)  is  the 
only  explicit  display  quantity,  the  (model)  display  vector  contains 
the  two  components: 


y^(t)  = e(t) 
y2(t)  = e(t) 


3*2.2  Performance  Requirements 


Performance  requirements  are  specified  in  terms  or  a "cost" 
that  is  to  be  minimized.  In  the  case  of  a single-variable 
tracking  task,  a good  match  between  model  and  measured  performance 
is  obtained  if  we  assume  that  the  pilot  is  attempting  to  minimize 
a weighted  sum  of  mean-squared  tracking  error  plus  mean-squared 
control  rate.  For  a system  perturbed  by  zero-mean  disturbances, 
the  cost  [J]  is  defined  as 
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where  G is  the  relative  "cost  weighting"  associated  with  mean- 
squared  control-ra-»-  activity.  (We  may  assign  a relative  weighting 
of  unity  to  mean-squared  error  with  no  loss  of  generality.)  One  of 
the  consequences  of  assigning  a cost  weighting  to  lontrol-rate, 
rather  than  to  control  displacement,  is  the  introduction  of  a 
first-order  lag  in  thi:  predicted  pilot  describing  function.  We 
refer  to  the  time  constant  of  this  lag  networ)  as  the  "motor  time 
constant" . 


For  a more  realistic  flight  task  in  which  the  pilot  mani- 
pulates a number  of  controls  to  regulate  a number  of  flight 
variables,  the  cost  will  be  defined  as 
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where  is  the  number  of  display  varia)>les;  the  number  of 
controls;  and  q^,  r^,  and  are  the  weighting  coefficients 
associated  with  the  various  display,  co:itrol,  and  control-rate 
variables,  respectively. 
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Reasonable  model  predictions  can  be  obtained  if  weighting 
coefficients  are  selected  on  the  basis  of  maximum  allowable 
deviations  (or  "limits”)  for  problem  variables  [3].  A unit 
"cost"  is  associated  with  a given  variable  when  the  magnitude 
of  the  "error"  (i.e.,  deviation  from  trim)  is  equal  to  the 
nominal  limit.  Thus,  the  weighting  coefficient  for  each 
variable  in  the  cost  equation  is  computed  simply  as  the  inverse 
of  the  square  of  the  corresponding  limit. 

3.2,3  Pilot  Limitations 

As  noted  above,  model  parameters  that  explicitly  relate  to 
pilot  limitations  are  effective  perceptual  time  delay,  visual  and 
motor  noise  processes,  and,  in  the  case  of  a single-variable  task, 
the  motor  time  constant.  Numerical  values  for  these  parameters 
have  been  found  to  be  remarkably  invariant  for  a variety  of 
tracking  situaf'  ns  [1,  2],  Time  delays  have  generally  fallen 
within  tiu-3  range  of  0.15  - 0.2  seconds,  and  motor  time  constants 
in  tlic  range  of  0.0  8 to  0.1  seconds  are  typical. 


A very  simple  model  accounts  for  pilot  remnant  in  situations 
where  control  gain  is  optimized  and  where  visual  threshold  and 
saturation  effects  are  negligible.  In  this  case,  both  observation 
and  motor  noise  processes  appear  to  scale  with  the  variance  of 
tlic  corresponding  problem  variables.  Thus,  in  the  single-variable 
tracking  task,  the  variance  of  the  injected  white  noise  processes 
are 
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where  P^#  P|,  and  P^  are  the  constants  of  proportionality  (i.e., 
"noise/signal  ratios”)  for  error  and  error~rat.'i  perception  and  tor 
motor  output,  respectively. 

Consistent  values  have  been  fjund  for  these  noise/signal 
ratios  for  a variety  of  labcratory  tasks.  Ratios  of  .Olir 
(i.e.,  -20  db)  are  typical  for  both  P^  and  P*,  and  a motor 
noise/signal  ratio  of  .003r  (-25  dD)  is  typical.  The  relative 
invariance  of  these  parameter*^  in  idealized  laboratory  situations 
suggests  that  the  noise/sig  -a.  ratios  reflect  fundamental  limita- 
tions on  human  capability  to  process  information. 

The  expression  for  observation  noise  variance  is  modified 
as  follows  to  account  for  additional  display-related  limitations 
typical  of  realistic  flight  tasks t 
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where  y^  denotes  the  i^j^  display  variable,  a^^  is  the  effective 
perceptual  threshold  associated  with  the  i^^^  variable,  K(a^,  a^^^) 
is  the  corresponding  describing-function  gain,  and  is  the 

variance  of  the  "residual  nois*;"  15].  The  effect  of  the 
threshold-related  gain  is  to  iiicrease  the  observation  noise 
without  limit  as  the  signal  deviation  becomes  small  with  respect 
to  the  "threshold",  thus  rendering  information  unobtainable  from 
the  display.  The  residual  noise  variance,  on  the  other  hand, 
represents  a constant  increment  of  observation  noise  resulting 
from  non-threshold  types  of  display  limitations. 
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Laboratory  studies  indicate  that  an  effective  threshold  of 
from  0.05  - 0.1  degrees  visual  arc  is  typical  for  perception  of 
indicator  displacement  and  a threshold  of  about  0,2  degrees/sec 
visual  arc  is  appropriate  for  velocity  perception  in  a tracking 
task  [4].*  If  the  pilot  is  required  to  maintain  a display 
indicator  at  some  non-zero  distance  from  a reference  indicator, 
the  residual  noise  variance  may  be  set  equal  to  the  square  of 
the  distance  (in  experimental  units)  from  the  indicator  to 
reference  [5]. 

3.2.4  Pilot  Optimization 

The  pilot  it>  assumed  to  adapt  his  control  strategy,  subject 
to  the  limitations  described  above,  to  minimize  the  cost.  The 
adaptive  portion  of  the  pilot's  strategy  is  represented  by  (a) 
an  optimal  predictor  that  partially  compensates  for  the  inherent 
time  delay,  (b)  an  optimal  estimator  to  obtain  the  best  estimate 
of  the  state  vector  x (t) , and  (c)  a set  of  optimal  gains  acting 
on  the  best  estimate  of  the  state  vector,  the  output  of  which 
is  a "commanded"  control  signa.l . This  output  is  modified  by  a 
first-order  low-pass  filter  (whose  time  constant  is  the  "motor 
time  constant"  mentioned  above)  to  reflect  inherent  limitations 
on  pilot  bandwidth. 


^Tbese  values  have  been  derived  mostly  from  studies  of  pilot 
remnant  in  single-variable  tracking  tasks  and  are  intended  to 
allow  a valid  statistical  linearization  of  apparent  threshold- 
related  phenomena  inherent  in  these  experimental  situations. 
The  reader  should  not  expect  these  "thresholds"  to  be 
necessarily  identical  to  other  measures  of  visual  threshold 
that  appear  in  the  literature. 
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3.2.5  Model  Outputs 


The  primary  output  varic.  I.3  of  this  pilot/vehicle  model  is 
the  pilot's  control  signal#  u(t).  Once  the  characteristics  of 
this  signal  have  been  determined#  a set  of  relatively  straight- 
forward computations  yields  the  pilot's  describing  function# 
remnant  spectrum#  and  various  measures  of  overall  system  per- 
formance (such  as  mean-squared  error  and  mean-squared  control 
activity) . 

3.3  Effects  of  Vibration 

Whole-body  vibration  will  generally  cause  unwanted  control 
inputs  due  to  biomechanical  coupling  between  the  vibration  source 
and  the  control  stick.  In  part#  then#  the  effect  of  vibration 
may  be  reflected  simply  as  an  additional  disturbance  to  the  vehicle 
as  diagrammed  in  Figure  8.  in  many  situations#  however,  the 
dominant  effect  of  vibration  is  to  interfere  with  the  sensory# 
central-processing,  and  motor  response  systems  employed  by  the 
pilot  in  performing  the  tracking  task.  Effects  of  this  sort  are 
modeled  as  changes  in  model  parameters  related  to  pi lot- limitations 
(i.e.,  time  delay#  motor  lag#  and  visual  and  motor  noise). 

The  treatment  of  various  sources  of  vibration-related 
performance  decrement  are  discussed  below. 

3.3.1  Vibration  Feedthrough 


Vibration-Correlated  Control  Inputs 

"Vibration  feedthrough"  is  the  term  ws  use  to  describe 
control  inputs  that  arise  from  (basically  linear)  biomechanical 
coupling  between  the  vibration  source  and  the  control  stick. 


Data  analysis  performed  in  this  study  (see  Appendix  B)  indicates 
that  the  impedance  model  diagrammed  in  Figure  9 may  be  useful 
for  predicting  the  interaction  of.  control-stick  dynamics  and 
feedthrough. 


We  define  the  transfer  impedance  ZT(jw)  as  the  force  trans- 
mitted to  a rigid  control  stick.  Thus, 


ZT 


ZS  * 


where  is  the  resulting  force,  is  the  platform  acceleration, 
and  ZR  is  the  mechanical  impedance  of  the  control  stick.  (For 
ease  of  notation,  wo  shall  delete  the  argument  (jw)  from  frequency- 
domain  models.)  The  output  impedance  ZO  is  defined  as  the  ratio 
of  the  rigid-stick  force  to  the  free-stick  displacement: 


ZO  J 
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"P  ' 

'zs  = « 

£v 
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ZS  = 0 

where  is  stick  disi^lacement . 


The  impedance  parameters  ZT  and  ZO  are  assumed  to  reflect 
biodynamic  properties  of  the  pilot  and  of  the  pilot-cockpit 
interfacf?  and  are  assumed  to  be  independent  of  control-stick 
characteristics  as  well  as  of  parameters  of  the  vibration  inputs. 
The  transfer  functions  for  vibration  feedthrough,  relating  stick 
force  and  stick  displacement  to  platform  vibration,  are; 


^v  _ ZT  • ZS 
- ZO  + ZR 


a ZO  + ZS 

P 


(2) 

(3) 
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vowin'':; 


ZT  « Transfer  Impedance 
70  = Output  Impedance 
7S  = Stick  Impedance 
= Platform  vibration 
0^  = Stick  Displacement 
= Stick  Force 


Figure  9.  Impedance  Model  for  Vibration  Feedthrough 
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In  order  to  obtain  numerical  predictions  of  vibration 
feedthrough,  v/e  must  hava  a means  of  deriving  the  feedthrough 
impedance  functions  ZT  and  ZO.  These  functions  will  in  general 
depend  on  the  structure  of  the  pilot's  seat,  the  interface 
between  the  pilot  and  the  seat,  the  interface  between  the  seat 
and  the  platform,  and  the  physical  characteristics  of  the  pilot. 

Efforts  are  currently  underway  to  develop  and  refine  models 
for  predicting  biodynamic  parameters  [6] , Existing  models  for 
vertical-axis  vibration  feedtiirough  range  from  simple  second-order 
systems  to  high-order  representations  containing  detailed  models 
of  the  neuromuscular  system  [7-8] . At  present,  however,  develop- 
ment of  biodynamic  models  is  not  sufficiently  advanced  to  allow 
reliable,  numerical  predictions  of  ZT  and  ZO.  Accordingly,  we 
must  rely  upon  measurements  obtained  in  the  specific  biodynamic 
environment  of  interest. 

Once  the  impedance  parameters  ZT  and  ZO  have  been  quantified 
and  the  vibration  spectrum  defined,  the  effects  of  vibration- 
correlated  control  inputs  on  tracking  performance  can  be  computed 
for  various  control-stick  characteristics  and  for  various  plant 
dynamics  as  described  below. 

Effect  on  Tracking  Performance 

Vibration  feedthrough  will  increase  tracking  error  because 
of  the  unwanted  inputs  applied  to  the  control  system.  Prediction 
of  this  error  increment  is  particularly  simple  if  the  vibration 
spectrum  contains  most  of  its  power  at  frequencies  beyond  the 
response  bandwidth  of  the  control  system.  In  this  case,  the 
increment  in  tracking  error  is  computed  as  the  open-loop  response 
of  the  control  system  to  the  feedthrough. 
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Assume  that  the  control  system  is  single-input/  single- 
output with  the  transfer  function  Y . and  let  K represent  the 

w W 

control  gain  in  terms  of  electrical  input  to  the  system  per  unit 
of  stick  displacement.  The  error  variance  predicted  by  feedthrough 
is  then 


0 


ZT 

J 

z6  V z "5^ 

•^cipap  dw 


(4) 


Since  the  vibration  input  is  assumed  to  be  uncorrelatad  with  the 
tracking  input,  this  error  term  can  be  added  to  the  error  variance 
arising  from  the  tracking  disturbance  to  yield  total  error 
variance.  This  expression  can  also  be  used  to  predict  various 
error  terms  of  interest  in  a multi-output  system,  where  the 
vehicle  transfer  function  is  now  represented  as  to  denote 
the  response  of  tlie  i^^  output  variable  to  the  control  input. 

If  the  vibration  spectrum  contains  significant  power  within 
the  response  bandwidth  of  the  pilot/vehicle  system,  the  pilot  will 
be  able  to  compensate  for  a portion  of  the  feedthrough-related 
tracking  error.  In  this  case,  feedthrough  error  would  be  less 
chan  that  predicted  by  the  open-loop  approximation  suggested 
above.  An  alternate  method  for  predicting  the  effects  of  feed- 
through would  bc'  to; 


(a)  Compute  the  spectrum  of  the  feedthrough-related 
stick  motion,  given  as 


Cv  ®V 


ZT 


ZO  + ZS 
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(b)  Approximate  this  spectrum  by  filtered  white  noise  and 
find  a suitable  state-variable  representation  for  this 
process , 


(c)  Treat  the  feedthrough  as  an  additional  tracking  input, 
and  use  the  state-variable  model  described  in  Section 
3.2  to  predict  response  behavior  simultaneously  to 
the  tracking  and  vibration  inputs. 


If  the  system  under  study  is  a moving-base  simulation  or  an 
actual  flight,  prediction  of  feedthrough  effects  is  further  compli- 
cated if  appreciable  vehicle  motions  (for  example,  high-frequency 
bending  modes)  are  induced  by  vibration-related  control  inputs. 
Such  motions  are  likely  to  be  fed  back  to  the  control  stick 
tJirough  biomechanical  coupling.  To  predict  feedthrough-related 
error  in  this  situation,  the  vehicle  dynamics  will  have  to  be 
augmented  to  account  for  that  additional  loop  closure. 


3.3.2  Visual  Effects 


Calaulatiort  of  Relative  b'ye-Dieplay  Motion 

Platform  vibration  will  in  general  produce  relative  motion 
between  tlie  head  and  the  platform  as  well  as  front-back  rotation 
of  the  head.  Translation  and  rotation  of  the  head  will  combine 
to  produce  relative  motion  between  the  head-point-of-regard  and 
the  display  which,  to  the  extent  that  the  oculomotor  system 
cannot  compensate,  will  produce  retinal  motion  of  the  display 
( i . e.  , blur ) . 


I 
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Once  the  vibration  environment  has  been  specified  and  the 
biodynamic  response  properties  of  the  head  determined,  rms  relative 
displacement  and  velocity  of  the  target  image  may  be  computed  as 
follows.  (Frequency-domain  representation  is  implied.) 

Let  = translational  acceleration  of  the  head  (g) 

: = translational  acceleration  of  the  display  (g) 

2 

olq  = rotational  acceleration  of  the  head  (rad/sec  ) ; 

Op  = platform  acceleration  (g) 
ji.  D = viewing  distance  (inches) 

■ ' 2 

f Acceleration  in  inch/sec  of  the  head-point-of-, regard , relative 

to  the  display,  is  then 

= 306.4  - a^)  + D • 

which  may  be  expressed  in  terms  cf  transfer  functions  as 


( 

[ where  the  constant  386.4  represents  the  acceleration  of  gravity, 

\ at  sea  level,  in  inches/sec^.  Typical  values  for  the  display/ 

j;  platform  transfer  function  are  unity  for  a display  rigidly  attached 

I to  the  accelerating  platform  (as  was  the  case  in  this  study) , and 

I zero  for  a completely  isolated  display  (as  is  sometimes  the  case 


Since  the  oculomotor  system  will  compensate  for  target 
motions  at  low  frequencies,  the  above  expression  must  be  modified 
by  a high-pass  filter  to  provide  an  estimate  of  the  relative 
acceleration  of  the  eye  point-of- regard.  Thus, 

“du  ' “dh  • ’'he 

Tlie  above  expressions  may  be  combined  to  yield  the  spectrum  of  the 
acceleration  of  the  image  motion: 


‘J'a^je^de 


386.4 


+ D 


Displacement  and  velocity  variances  are  then  given  by  the  following 
integrations : 


„ 2 _ r “de°^de 

‘J — 7“ 


dw 


(5) 


2 _ r *^de<^de 

^de  “i  7“ 


doj 


(fi) 


As  we  do  not  presently  have  reliable  theoretical  modeli^  for 
relating  head  translation  and  rotation  to  platform  acceleration, 
the  quantities  a^/otp  and  oi0/ap  mus':  be  determined  experimentally 
in  biodynamic  environments  of  interest.  The  transfer  function 
a^/Op  is  determined  by  the  interface  between  the  display  and  the 
vibrating  platform. 
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A simple  model  can  be  suggested  for  the  transfer  function 
*^he*  Laboratory  studies  have  shown  that  the  eye  can  follow  target 
motion  that  contains  frequencies  below  2 Hz.  Response  characteris- 
tics show  a very  rapid  falloff  between  2 and  4 Hz,  and  the  tracking 
effectiveness  of  the  eye  at  frequencies  above  4 Hz  is  essentially 
nil  [0-91,  Since  the  parameter  is  the  complement  of  the 
transfer  function  describing  oculomotor  response,  may  be 
approximated  reasonably  well  by  a step  function  in  the  frequency 
domain  that  has  zero  gain  from  d-c  to  about  4 Hz  and  unity  gain 
above  4 Hz. 

Effect  on  Tracking  Performance 

Appropriate  modifications  to  the  parameters  of  the  pilot/ 
vehicle  model  of  Section  3,2  must  be  made  in  order  to  account  for 
the  effects  of  relative  image  motion  on  tracking  performance.  It 
is  not  clear  from  theoretical  considerations,  however,  just  what 
these  modifications  should  be.  If  vibration-induced  motion 
affected  only  the  error  indicator  and  not  the  reference  indinator, 
we  could  treat  image  motion  as  a noise  added  to  the  error  display; 
in  this  case,  the  mean-squared  image  displacement  and  velocity 
would  be  treated  as  residual  noise  variances  and  would  be  included 
in  the  expression  for  observation  noise  as  indicated  in  eauation 
(1)  . 


Since  vibration-induced  blur  causes  simultaneous  retinal 
motion  of  both  the  error  indicator  and  the  corresponding  reference 
indicator,  it  is  not  o)-»vious  how  blur  degrades  the  pilot's  ability 
to  distinguish  the  position  and  velocity  of  the  error  indicator 
relative  to  the  reference  indicator.  We  must  therefore  rely  on 
experimental  results  to  serve  as  a guide  for  modeling  the  effects 
of  blur. 


30 


The  effect  of  vibration-rexated  blur  appears  to  depend  on 
the  nature  of  the  vibration  spectrum.  In  a study  conducted  by 
Allen  et  al.  usina  sinusoidal  vibration^  the  test  subjects 
reported  that  their  ability  to  allign  error  and  reference  lines 
seemed  to  be  unimpaired  by  blurring  [7].  Allen  et  al.  commented 
that  the  luminous  density  of  the  blurred  images  of  the  error  and 
reference  lines  were  greatest,  near  the  extreme  of  travel,  thereby 
providing  sharp  edges  for  allignment, 

Since  a auasi-random  (sum  of  five  sinusoids)  vibration  was 
used  in  the  current  study,  the  blurred  images  were  expected  to  be 
somewhat  more  difficult  to  allign.  Therefore,  the  effect  of  blur 
was  approximated  by  setting  the  residual  noise  variance  associated 
with  error  displacement  equal  to  the  displacement  variance  calculated 
);'/  oouation  (B)  . Model  predictions  obtained  with  this  t.reatment 
of  blur  were  consistent  wit)',  experimental  results.  On  the  otlier 
hand,  setting  residual  noise  on  rate  perception  equal  to  image 
velocity  as  calculated  by  equation  (f>)  yielded  a predicted  pilot 
remnant  spectrum  that  differed  appreciably  from  that  observed 
oxnerimontally  . A better  model  match  was  obtained  vrith  image  rate 
ignored;  i.e.,  with  rate  residual  noise  set  to  zero. 

Until  we  obtain  a better  theoretical  understanding  of  how 
vibration-induced  image  motion  influences  perceptual  processes 
used  in  tracl^ing,  we  suggest  that  the  effects  of  blur  be  accounted 
for  by  a residual  noise  variance  associated  with  perception  of 
error  displacement.  This  variance  is  to  be  ecual  to  the  variance 
of  tliG  eye  point-of-reqard  displacement  relative  to  tlie  display. 
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3.3.3  Motor  Effects 


In  addition  to  inducing  control  motions  correlated  with  the 
vibration  input  as  discussed  above,  whole-body  vibration  can 
interfere  with  the  pilot's  ability  to  execute  the  desired  response 
to  the  tracking  input.  Interference  of  this  sert  has  been  thought 
to  arise  largely  from  degradation  of  sensory  f jedback  information 
used  in  the  regulation  of  control  activity  [7,  8,  10] , leading  to 
measureable  effects  such  as  increased  motor-related  pilot  remnant 
and  a more  sluggish  tracking  response.  Given  the  existing  structure 
of  the  state  variable  pilot/vehicle  model,  one  would  expect  these 
motor  effects  to  be  mo^-ided  by  one  or  more  of  the  follov/ino  para- 
meter modifications s (a)  increased  motor  noise,  (b)  increased 
time  delay,  and  (c)  increased  motor  time  constant. 


For  the  biodynamic  environment  explored  in  this  study, 
motor-related  effects  of  vibration  can  be  accounted  for  by  increases 
in  the  motor  noise/signal  ratio  and  in  time  delay;  no  change  is 
required  for  the  motor  time  constant.  (See  Appendix  D for 
detailed  discussion  of  the  modeling  performed  in  this  study.) 

We  suspect  that  the  increased  motor  noise  reflects  an  increase 
in  the  noise  associated  with  proprioceptiye  and  kinematic  sensory 
pathways,  and  that  the  pilot's  adaption  to  this  additional  noise 
results  in  describing-function  changes  that  reauire  an  increased 
time  delay  for  adequate  model-matching.*  Further  research  is 
needed  to  determine  (a)  how  the  injected  motor  noise  varies  with 
vibration-related  quantities  in  a vibration  environment,  and  (b) 
how  the  effective  time  delay  relates  to  vibration  and/or  noise 
parameters . 

*The  re"ader  is  remTnded  that  "time  delay"  refers  to  a model  para- 
meter that  may  reflect  the  effects  of  a number  of  psychophysiological 
phenomena.  Increased  time  delay  is  reguired  also  by  existing 
frequency-domain  models  to  adequately  mimic  the  effects  of  vibra- 
tion [ 8 ] . 
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Analysis  recently  performed  on  manual  tracking  data 
obtained  in  a lateral  vibration  environment  suqqested  that 
motor  noise  could  be  considered  to  scale  with  overall  control 
activity;  i.o.,  vibration  feedthrough  plus  response  to  tracking 
error  [10].  Analysis  of  longitudinal-axis  vibration  performed 
in  the  current  study  does  not  support  this  hypothesis,  however. 

An  alternative  hypothesis  more  consistent  with  the  existing  data 
base  is  that  the  motor  noiso/siqnal  ratio  varies  v;ith  lim)-'  (o.g., 
slioulder  and/or  elbow)  motion. 

Studies  of  trackino  in  a vibration  environment  --  including 
the  study  documented  in  this  report  — indicate  tliat  the  following 
spring-mass-damper  model  is  sufficient  to  relate  shoulder  motions 
to  platform  acceleration  [7] ; 


l£(s)  - 

K + Bj,  + M s2 

_ 1 H-  (B/K)S 

r + rn7KTs+  (m/k  ) s2‘ 


where  K,  B,  and  M represent  spring  constant,  damping  factor  and 
mass,  respectively.  Once  the  ratios  B/K  and  M/K  have  been  deter- 
ruined  for  pilots  of  a particular  anthropometric  structure  situated 
in  a 5'.pecific  biodynamic  environment,  the  variance  of  the  shoulder 
acceleration  (or  of  any  linear  transformation  of  acceleration)  can 
be  predicted  for  various  vibration  spectra  as  follows; 
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where  o is  the  variance  of  t).f,  quantity  of  interest  (e.g.  , the 
first  derivative  of  the  acceleration)  and  F is  the  transformation 
process  (e.g.,  differentiator)  expressed  as  a transfer  function. 
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In  summary,  current  understanding  of  the  effects  of  vibration 
on  tracking  behavior  is  not  sufficient  to  allow  us  to  specify  ana- 
lytic relationships  between  vibration  parameters  and  changes  in 
motor-related  model  parameters;  we  must  therefore  detenr.ine  model- 
parameter  changes  from  a "calibration  experiment"  performed  in  a 
specific  biodynamic  environment.  Hopefully,  the  results  of  such 
an  experiment  will  be  valid  for  a variety  of  vibration  environments 
and  a variety  of  control-system  response  characteristics, 

3.3.4  Central-Processing  Effects 


Long-term  exposure  to  vibration  can  be  expected  to  lead  to 
stress  effects  such  as  fatigue,  discomfort,  anxiety,  etc.,  which 
would  further  degrade  the  pilot's  ability  to  process  tracking 
information.  Substantial  research  is  needed  to  develop  models 
for  long-term  effects,  as  studies  of  manual  tracking  in  vibra- 
tion environments  have  dealt  primarily  with  tracking  tasks  of 
a few  minutes  duration.  Consequently,  wc:  shall  not  attempt  to 
account  for  these  effects  in  the  design  guide  presented  in  the 
following  section. 
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4.  DESIGN  GUIDE 


Guidelines  for  the  design  of  control  sticks  in  a vibration 
environment  can  take  a number  of  forms,  depending  on  the  intended 
range  of  applicability.  The  more  restricted  the  range,  the  easier 
it  is  to  present  the  design  guide  in  a form  that  is  immediately 
unseable.  For  example,  if  we  are  concerned  with  the  specification 
of  a single  stick  parameter  (say,  control  gain)  for  a specific 
tracking/vibration/biodynamic  situation,  the  guide  can  consist  of 
a single  curve  of  performance  (e.g.,  mean-squared  tracking  error) 
versus  the  parameter  of  interest.  On  the  other  hand,  if  the  guide 
is  to  apply  to  a range  of  tracking  and  biodynamic  situations,  it 
will  ideally  take  the  form  of  an  analytic  model  (or  set  of  models) 
having  parameters  that  are  either  independent  of  the  specifics  of 
the  situation  or  that  can  be  readily  determined  on  a theoretical 
basis.  In  this  section  we  shall  consider  the  proposed  design  guide 
from  both  the  general  and  restricted  views. 

The  proposed  design  guide  is  based  on  the  model  structure 
outlined  in  the  preceding  section  of  this  report.  In  the  most 
general  sense,  the  model  itself  ie  the  guide;  that  is,  the  model 
is  to  be  applied  in  a specific  situation  cf  interest  to  explore 
the  consequences  of  one  or  more  candidate  control-sticlc  confi- 
gurations . 

The  status  of  biodynamic  modeling  prohibits  a fully  general 
design  guide  at  this  time.  In  many  cases,  analytic  forms  for 
biodynamic  response  mechanisms  cannot  reliably  be  presented. 
Furthermore,  where  such  forms  exist,  we  do  not  know  how  to  quantify 
parameters  without  the  benefit  of  experimental  data  obtained  in  the 
specific  biodynamic/vibration  environment.  In  addition,  we  do  not 
know  the  functional  relationship  between  vibration  parameters  and 
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changes  in  parameters  of  the  pllot/vehlcle  model.  Consequently, 
the  design  guide  discussed  here  must  be  restricted  to  the  bio- 
dynamic and  vibration  environments  explored  in  this  study  as 
described  in  Section  2.  Nevertheless,  the  highly  analytic 
pilot/vehicle  model  Incorporated  in  the  design  guide  should 
allow  one  to  explore  tracking  tasks  and  control  configurations 
that  have  not  been  studied  experimentally. 

Quantification  of  model  parameters  for  the  design  guide  is 
discussed  belo’*'  in  Section  4.1.  Numerical  examples  are  presented 
in  Section  4.2  to  validate  the  model  and  to  provide  simple 
quantitative  guidelines  for  the  selection  of  specific  control 
parameters  in  specific  situations. 

4.1  Quantification  of  Model  Parameters 

Model  parameters  are  categorized  as  follows  for  convenience 
of  discussion:  (1)  data-dependent  parameters,  (2)  problem 

variables,  and  (3)  remaining  pilot  parameters. 

For  the  most  part,  parameters  of  the  design  guide  that 
relate  directly  to  the  vibration  and  biodynamic  environments  are 
empirical  parameters  that  have  been  quantified  on  the  basis  of  the 
experimental  study  described  in  Section  2.  These  parameters  are 
considered  to  be  data-dependent;  that  is,  we  have  no  theoretical 
basis  for  reliably  predicting  how  parameter  values  will  change 
for  different  biodynamic/vibration  environments.  On  the  other 
hand,  parameters  that  describe  the  tracking  task  and  the  control 
configuration  can  be  considered  as  variables  of  the  design  guide. 
The  remaining  pilot-related  model  parameters  are  either  relatively 
Invariant  or  can  be  determined  on  a theoretical  basis. 
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Data  dependent  parameters,  problem  variables,  and  remaining 
pilot  parameters  are  discussed  separately  below.  The  reader  is 
assumed  to  be  familiar  with  the  model  structure  and  parameter 
definitions  presented  in  Section  3 of  this  report. 

4.1.1  Data-Dependent  Parameters 

Vibration  Environment 

This  design  guide  applies  to  a z-axis  vibration  disturbance 
of  about  .26  g rms  consisting  of  the  sum  of  five  sinusoids.  The 
average  vibration  spectrum  measured  in  this  experimental  study  is 
shown  below  in  Table  2.* 


Table  2 

Experimental  Vibration  Spectrum 

Frequency  Power 

(Hz)  (rad/sec)  (dB) 


2 

12.5 

-17.7 

3 

21.0 

-18.1 

5 

31.7 

-19.2 

7 

43.8 

-19.7 

10 

63.  3 

-20.1 

Feedthrough 

Impedance  Parameters 

Stick  displacement  is  assumed  to  be  correlated  with  platform 
vibration  as  follows; 

*See  Appendix  A ior  a definition  of  spectral  quantities. 
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ZT 

ZO  + ZS 


(8) 


where  the  transfer  Impedance  ZT  and  the  output  Impedance  ZO  are 
as  shown  In  Figure  10  for  center  and  side  stick  locations.  An 
amplitude-ratio  of  0 dB  signifies  1 lb  of  stick  force  per  g of 
platform  acceleration  for  ZT  and  1 lb  of  stick  force  per  inch 
of  stick  displacement  for  ZO. 


Relative  Eye-Dieplay  Motion 


The  spectrum  of  relative  motion  between  the  eye-point-of- 
regard  and  the  display  may  be  estimated  as  follows: 


o^de^de 


1386.4 


!h . D . ^ 

“p  “p/  “p 


he 


^apOp  (c,) 


where  the  transfer  functions  relating  vertical  head  acceleration 

to  platform  acceleration  (uj^/Up)  and  head  rotation  to  platform 

acceleration  (a„/a  ) are  as  shown  in  Figure  11.  The  expression 
up  2 

) is  dimensionless,  whereas  a„/a_  has  units  of  rad/sec  )/g. 

« P 

The  filter  function  is  assumed  to  have  zero  transmission  from 
d.c.  to  4 Hz  (25  rad/sec)  and  unity  transmission  for  frequencies 
above  4 Hz. 


Viewing  distance  D and  display/platform  transfer  may  be 
treated  as  problem  variables,  if  desired.  For  the  data  analysis 
performed  in  this  study,  D was  assumed  to  be  30  inches,  and  the 
transfer  function  (oij/a  ) was  unity  to  reflect  the  fact  that  the 
display  was  rigidly  attached  to  the  platform. 
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Once  the  spectrum  of  eye-display  motion  is  computed,  the 
integration  defined  by  equation  (3-5)  may  be  performed  to  yield 
the  variance  of  the  relative  motion  between  the  display  and  the 
eye-point-of-regard,  Thn  residual  noise  variance  associated 
with  perception  of  tracking  error  (a  model  parameter)  is  set 
equal  to  this  derived  measure. 

Vibration-Related  Pilot  Parametera 

A time  delay  of  0.2  seconds  and  a motor  noise/signal  ratio 
of  0.04  (relating  the  variance  of  the  injected  white  motor  noise 
to  the  variance  of  the  predicted  control  force)  are  assigned  on 
the  basis  of  measurements  obtained  in  this  study.  These  parametf  ■ 
values  are  suggested  only  for  the  vibration  and  biodynamic  environ- 
ments described  in  this  report.  (Values  of  0.15  seconds  and  0.004 
are  suggested  for  static  environments.) 


4,1.2  Problem  Variables 


Task  Parametera 

Previous  application  of  the  pilot/vehicle  model  has  shown 
that  a good  match  between  model  outputs  and  experimental  measure- 
ments can  be  obtained  with  a consistent  set  of  pilot-related 
parameters  for  a variety  oi  laboratory  tracking  tasks  [1,  2], 

On  this  basis,  we  shall  assume  that  the  design  guide  described 
here  can  be  expected  to  give  reasonably  accurate  predictions  for 
tracking  tasks  usxng  vehicle  dynamics  and  input  spectra  other 
than  those  explored  in  this  study.  Application  is  restricted, 
however,  to  the  static  situation  and  to  the  vibration/biodynamic 
environment  explored  in  this  study. 


47 


II  IN-  IHl1l“' 


n itniit  • riif  I'i 


Dieplay •Related  Parameters 

The  design  guide  should  be  valid  for  a variety  of  display 
designs.  Display  dynamics  are  included  in  the  description  of 
system  dynamics  as  discussed  in  Section  3.2.1.  Limitations  of 
visual  resolution  may  also  be  treated  by  appropriate  selection 
of  effective  "thresholds".  Typically,  position  and  rate  informa- 
tion will  be  obtained  from  each  moving  display  element.  On  the 


basis  of  previous  studies  of  pilot  remnant  [4] , we  suggest  that 
threshold  levels  be  computed  as  follows: 


D 0. 05 

s * vrn 


D 0*2 

G * TT7S 


where  a.,  and  a*  are  the  thresholds  associated  with  displacement 

y y 

and  rate,  expressed  in  the  desired  experimental  units  (e.g., 
volts);  D is  the  distance  to  the  display  (e.g.,  inches),  and  G 
is  the  display  gain  expressed  as  units  of  display  deflection  per 
experimental  unit  (e.g^ , inches/volt) . 

Application  of  the  design  guide  is  suggested  for  situations 
in  which  the  display  format  is  symbolic,  rather  than  pictorial, 
and  in  which  saturation  effects  are  not  important. 

Control’-Ctiok  Parameters 


The  design  guide  is  intended  especially  to  explore  the 
effects  of  control-stick  parameters  on  performance  in  the  vibration/ 
biodynamic  environment  described  in  Section  2.  Application  to 
stick  configurations  not  explored  in  this  study  is  demonstrated 
in  Section  4.2. 


Although  any  linear  relationship  between  stick  force  and 
stick  displacement  may  be  considered,  we  suspect  that  the  following 
second-order  impedance  functions  will  be  adequate  in  most  cases: 


“‘S'®’  “ 

where  ZS  is  the  impedance  relating  force  in  pounds  to  displacement 

in  inches;  K is  the  spring  gradient  in  Ibs/inch;  B is  the 
s s 

damping  in  lbs/ (in/sec) ; and  M is  the  effective  mass  of  the 

s 

stick  in  pounds.  Stick  dynamics  are  included  in  the  description 
of  system  dynamics  as  demonstrated  in  Section  4.2. 


4.1.3  Remaining  Pilot  Parameters 

For  a well-trained  test  subject,  the  observation  noise/signal 
ratio  appears  to  reflect  the  level  of  attention  to  the  tracking 
task.  For  single-variable  tracking  tasks,  this  parameter  appears 
for  the  most  part  to  be  independent  of  the  specifics  of  the  tracking 
tasks  and  is  unaffected  by  the  presence  or  absence  of  vibration 
(at  least  for  the  short  exposure  intervals  that  have  been  studied 
experimentally).  A value  of  about  -20  dB  (i.e.,  O.Olir)  is  typical, 
although  a value  of  -21  dB  was  found  to  provide  a better  match  to 
the  data  explored  in  this  study. 

The  observation  noise/signal  ratio  will  generally  be  increased 
if  the  pilot  has  to  share  attention  among  a number  of  tasks  and/or 
among  a number  of  display  variables  within  a given  task.  A reasonable 
model  for  the  effect  of  attention  is  to  let  the  noise/signal  ratio 
vary  as  follows; 


'’vi  - 
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where  Pyj^  is  the  n -ise/signal  ratio  associated  with  the  i^j^  display 
variable,  is  the  fraction  of  attention  devoted  to  that  variable, 
and  is  the  nominal  value  associated  with  single-variable  tracking 
(say,  -20  dB)  (11]  . 

In  most  single-variable  laboratory  situations  that  we  have 
studied,  the  "cost"  that  the  pilot  is  assumed  to  minimize  •'-pears 
to  represent  a limitation  on  pilot  bandwidth  rather  than 
objective  performance  requirement.  Typically,  the  expression 
for  the  cost  is  a weighted  sum  of  mean-squared  tracking  error 
and  mean-squared  control-rate,  wh  re  the  relative  weighting  is 
chosen  to  provide  an  effective  "mt  tor  time  constant"  of  about 
0.1  second. 

If  excessive  control  forces  or  displacements  are  required, 
an  additional  weighting  coefficient  on  control  force  (or  displace- 
ment) is  needed.  The  f:>llowing  expression  for  cost  was  found  to 
match  the  data  (both  s ic  and  vibration  conditions)  analyzed 
in  this  study: 

J = 1.0  a\  + .03  aj  + g (13) 

2 2 2 
where  is  the  mean-squared  error  in  volts  , o is  the  mean- 

e ^2  u 

squared  force  in  pounds'^,  and  (g)  is  selected  to  provide  a motor 

time  constant  of  0.1  seconds.  Since  the  display  gain  used  in 

this  study  was  approximately  0.19  inches/volt,  the  relative 

weightings  assigned  to  error  and  to  control  force  imply  that  the 

test  subjects  considered  the  presence  of  a 1-inch  tracking  error 

as  "costly"  as  the  requirement  to  generate  about  30  pounds  of 

control  force. 
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As  discussed  in  Section  3.2.2,  a more  elaborate  cost 
expression  will  be  required  to  adequately  model  performance  in 
a realistic  multi-variable  flight  situation.  In  this  case, 
weighting  coefficients  reflect  objective  performance  requirements 
and  hardware  limitations  as  well  as  pilot  preference  [3]. 

Suggested  values  for  pilot-related  parameters  of  the 
tracking  model  are  summarized  in  Table  3.  The  absence  of 
numeric  entries  indicates  that  parameter  values  are  to  be 
computed  as  described  above. 

4.2  Numerical  Examples 

Three  examples  are  presented  to  demonstrate  application 
of  the  design  procedure  and  to  serve  as  guidelines  for  possible 
additional  experimentation  to  be  performed  by  AMRL  with  the  same 
set  of  control  devices  in  a similar  laboratory  setting.  The 
primary  control  variables  explored  in  the  following  examples 
are,  in  order,  (1)  spring  gradient,  (2)  electrical  stick  gain, 
and  (3)  stick  damping. 


4.2.1 


For  our  first  illustrative  example  of  the  design  procedure 
we  compare  model  predictions  with  the  results  of  the  experimental 
study  described  in  Section  2 of  this  report.  The  reader  will 
recall  that  the  spring  gradient  was  one  of  the  major  experimental 
variables,  along  with  stick  location  and  presence  or  absence  of 
vibration.  Because  of  equipment  constraints,  however,  stick 
damping  and  electrical  gain  also  had  to  be  varied  as  spring 
gradient  war  changed,  since,  in  general,  performance  will  be 
affected  by  a change  in  any  one  of  these  stick  parameters,  this 


Table  3 


Su^ested  Values  for  Pilot-Related  Model  Parameters 


Parameter 


Parameter  Value 
Static  I Vibration'^ 


2 2 
a~  + g af 

u ^ u 


0.20  sec 
0.04 


*Valid  only  for  the  vibration/biodynamic  environment  described 
in  Section  2. 

**Tracking  error  in  inches/  force  in  pounds,  and  g selected  to 
provide  a motor  time  constant  of  0.1  seconds. 

***Value  calculated  as  described  in  text. 
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example  will  not  lead  to  a simple  curve  of  performance  versus 
spring  gradient.  We  begin  with  this  example,  nevertheless,  as 
it  affords  the  best  means  for  quantitatively  validating  the 
design  procedure. 

The  relevant  control-stick  parameters  are  shov/n  in 
Table  ^ Spring  gradient  and  electrical  gains  were  determined 
in  this  experimental  study;  stick  damping,  stick  mass,  and  effec- 
tive hand  mass  were  obtained  from  documentation  of  the  preceding 
study  performed  at  AMRL  [7].  Vehicle  transfer  functions  relating 
system  error  (E)  to  disturbance  input  (I)  and  control  input  (C) 
are  reviewed  below,  with  all  variables  in  volts: 


f(S,  = 1 

(14a) 

- % 

(14h) 

Since  the  lowest  vibration  frequency  (2  Hz,  or  12.5  rad/sec) 
was  beyond  the  gain-crossover  frequency  of  the  pilot/vehicle 
system  (less  than  5 rad/sec) , tracking  response  and  vibration 
feedthrough  were  analyzed  separately  and  the  results  combined  to 
yield  predictions  of  system  performance  in  the  vibration 
environment. 

Response  to  Tvaaking  Input 

The  force-feel  properties  of  the  pilot/stick  interface  must 
be  accounted  for  in  order  to  model  adequately  the  effects  of 
control-stick  parameters  on  tracking  performance  in  static  and 
vibration  environments.  This  interaction  was  modeled  by  the 
following  impedance  function  relating  stick  output  C (in  volts) 
to  force  input  F (in  lbs) : 
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Table  4 


Control  Stick  Parameters 


Parameter 


Control  Stick 


Spring 

Light  Med  Hea^ 


Stiff 

Light  Med  Heav\ 


Electrical  gain, 
volts/lb 

1.2 

*^e 

Electrical  gain, 
volts/inch 

2.5 

Spring  gradient, 
Ibs/inch 

1 

2 

1 

Stick  damping,  [ 

lbs/ (inch/sec) 

O 

to 

Stick  mass,  lbs  j 

: 0.7 

Mass  of  hand  end  of; 
arm,  lbs  ! 

2.7 

1 

+ M^,  lbs  ; 

' 3.4 

027  .027 


(15) 


pis) 


K 

e 

2 


K. 


where  the  parameters  of  the  impedance  and  their  units  are  as 
expressed  in  Table  4.  Note  that  the  force  F is  not  necessarily 
the  force  applied  to  the  control  stick,  but  rather  the  force 
needed  to  drive  the  combination  of  the  stick  plus  the  inertia  of 
the  forearm. 

The  vehicle  transfer  functions  and  the  force-feel  impedance 
were  combined  to  yield  the  effective  control  system  shown  in  the 
flow  diagram  of  Figure  12.  The  reader  is  referred  to  Appendix  D 
for  a state-variable  representation  of  this  system. 

Since  the  model  analysis  described  in  this  application  of 
the  design  guide  was  performed  with  the  benefit  of  having  experi- 
mental data  in  hand,  pilot-related  model  parameters  were  selected 
to  provide  a good  match  overall  to  the  various  experimental  condi- 
tions. To  this  extent,  model  results  given  here  are  not  true 
"predictions".  Since  model  parameters  were  kept  fixed  for  all 
six  control-stick  configurations,  however,  the  diffevenoee  in 
performance  measures  that  the  model  exhibits  for  the  different 
configurations  represent  actual  predictions.  Furthermore,  model 
parameters  not  directly  affected  by  the  presence  of  vibration 
have  numerical  values  quite  similar  to  those  found  in  previous 
studies  of  tracking  in  non- vibration  environments. 

Parameters  relating  to  cost  functional,  time  delay,  and 
motor  noise  were  as  indicated  in  Table  3 for  the  vibration 
environment.  A somewhat  lower  observation  noise/signal  ratio 
(-21  dB)  than  that  shown  in  Table  3 was  used,  as  it  provided 
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I = Trackino  Disturbance,  volts 

F = Trackincr  Frror,  volts 

C = Control  Input,  volts 

F = Applied  Force,  pounds 

K = Flectrical  Ftick  Gain,  volts/inch 
e 

K^=  Sprinq  Gradient,  pounds/inch 

B,=  rtick  Dampinq,  pounds/ { in/sec) 
s 

Combined  Mass  of  Stick  and  Forearm,  pounds 


Fiqure  12,  Linear  Flow  Diaqram  of  Control  System 
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a better  match  to  both  static  and  vibration  data  obtained  in  this 
study,  "Threshold"  levels  of  approximately  0.14  and  2.7  volts 
were  computed  from  equation  (4-3)  for  displacement  and  rate 
perception,  respectively,  A display  gain  of  0,1875  inches/volt 
and  viewing  distance  of  30  inches  were  used  for  this  calculation. 
Computations  equivalent  to  those  implied  by  equation  (4-2)  led 
to  an  rms  residual  noise  on  error  displacement  of  0.44  volts.* 

Once  the  system  was  described  analytically  and  pilot  para- 
meters quantified/  the  pilot/vehicle  model  was  then  able  to 
predict  measures  of  system  performance  for  the  various  control 
configurations  in  the  absence  of  vibration  feedthrough.  Errors 
and  control  motions  due  to  feedthrough  were  computed  separately, 
as  discussed  below. 

Vibration  Feedthrough 

Vibration-correlated  error  and  control  scores  were  computed 
separately  for  the  side  and  center  stick  locations,  with  the 
results  then  averaged  together.  With  the  vehicle  dynamics  given 
as  in  equation  (14h)  vibration-correlated  error  was  predicted  as 


0 


2 

®v 


5 

I 

i=l 


ZT  16  „2 

Zd  + “Z^'S  “7  e 

U) 


(16) 


*See  Appendix  B for  detailed  discussion  of  the  computation  of 
relative  motion  between  eye-point-of-regard  and  display. 
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where  4>a  a*  represents  the  power  (in  g^)  contained  in  the  i 
sinusoidal  component  of  the  vibration  input  (Table  2)  t 
electrical  gain  of  the  stick  (Table  4)  , ZT  and  ZO  the  £eed~ 
through  impedance  functions  defined  in  Figure  10,  and  ZS  the 
stick  impedance.  The  stick  impedance  was  defined  as 


^s  2 

- 37S7J  ® * ® 


s ® * ■'s 


with  parameters 


and  as  defined  in  Table 
s 


4.  J'n 

expression  identical  to  equation  (16),  but  minus  the  quantities 
related  to  vehicle  dynamics,  was  used  to  predict  vibration- 
correlated  control  scores. 


I 


Performance  Prediotiona 


Predictions  of  tracking-related  and  vibration-correlated 
error  and  control  scores  are  shown  separately  and  combined  in 
Table  5.  The  feedthrough  component  of  tracking  error  is  small, 
accounting  for  less  than  five  percent  of  the  predicted  total 
error  in  all  cases.  Feedthrough  is  a much  more  prominent  aspect 
of  the  control  score,  however.  For  two  of  the  control-stick 
configurations  - light  and  medium  stiff  •-  vibration  feedthrough 
accounts  for  about  half  the  predicted  mean-squared  control 
displacement. 

Predicted  rms  error  and  control  scores  are  compared  with 
measured  performance  scores  in  Figure  13.  Each  experimental  data 
point  in  this  figure  represents  the  average  performance  of 
seven  test  subjects,  with  the  brackets  indicating  + 1 standard 
deviation  due  to  inter-subject  variability. 
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Control-Stick  Configuration 
Spring  Stiff 


Variable 


2 

°et 

(volts^) 

.250 

.229 

2 

0 0 

(volts^) 

.002 

.001 

2 

‘^etotal 

(volts^) 

.252 

.230 

°<^total 

(volts) 

.502 

.480 

2 

°ct 

(volts^) 

.231 

.273 

2 

®Cv 

(volts^) 

.045 

.016 

2 

°ctotal 

(volts^) 

.276 

1 

.289 

^ctotal 

(volts) 

.5  25 

.538 

2 

1 ^ 

2 i 

(pounds  ) 

.102 

1.19 

Light 

Stiff 

Med 

.227 

.229 

.009 

.010 

.236 

.239 

.485 

.489 

.473 

.541 

.448 

.434 

.921 

.975 

.960 

.988 

.228 

.231 

The  model  predicts  essentially  no  change  in  rms  error  score 
across  the  control-stick  configurations.  The  experimental  scores 
exhibit  some  variability  (which  is  statistically  significant) , 
but  no  consistent  trend.  Since  model  predictions  fall  within 
one  standard  deviation  of  the  corresponding  experimental  measure- 
ments, we  consider  the  agreement  between  model  and  measurement  to 
be  adequate  for  the  rms  error  score. 

Ag;.eement  is  much  closer  between  "predicted"  and  measured 
rms  control  scores,  however,  and  stiff-stick  scores  are  shown  to 
be  more  than  dov^ble,  on  the  average,  than  the  spring-stick 
ricores.  This  result  implies  the  spring  sticks  would  be 
the  better  choice  if  the  vehicle  exhibited  resonances  at  frequen- 
cies where  vibration  power  was  significant,  given  the  control- 
stick  configurations  explored  in  this  particular  study. 

The  relative  insensitivity  of  predicted  rms  error  with 
control-stick  configuration  shown  in  Figure  13  is  somewhat 
misleading  in  that  the  effects  of  vibration  have,  to  some  extent, 
counterbalanced  the  basic  performance  differences  across  control 
sticks.  As  shown  in  Appendix  C,  model  predictions  as  well  as 
experimental  measurements  show  improved  performance  with  the 
stiff  sticks  in  the  absence  of  vibration.  Vibration  causes  a 
greater  percentage  increase  in  rms  error  with  the  stiff  sticks, 
however,  which  tends  to  equalize  error  scores  across  spring 
and  stiff  sticks  when  vibration  is  present. 

4.2,2  Optimization  of  Electrical  Gain 


In  this  example  we  demonstrate  application  of  the  design 
method  to  the  optimization  of  a single  control-stick  parameter  in 
a particular  tracking  situation.  Specifically,  we  consider  the 


tracking  task  and  biodynainic/vibration  environment  defined  in  the 
preceding  example,  and  we  optimize  the  electrical  gain  of  the 
"light-stiff"  stick  as  defined  in  Table  4.  Because  we  are 
keeping  constant  all  remaining  parameters  of  the  task,  we  can 
derive  a design  curve  that  indicates  the  range  of  control  gains 
that  may  be  considered  optimal. 

As  in  the  preceding  example,  vibration  feedthrough  and 
tracking  response  are  computed  separately  and  combined.  From 
equation  (16)  we  note  that  the  variance  of  the  vibration- 
correlated  tracking  error  is  proportional  to  the  square  of 
the  stick  gain  when  all  other  system  parameters  are  held  fixed. 

We  may  thus  use  the  results  of  Table  5 to  predict  vibration- 
correlated  error  scores  (as  well  as  control  scores,  if  desired) 
to  any  control  gain  as  follows: 
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where  c 
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®o 


is 


the  error  score  computed  previously  for  the  gain  K^. 


Since  we  are  considering  a "stiff"  stick  in  which  stick 

force,  rather  tha^  ?tick  displacement,  is  considered  the  control 

quantity  of  interest,  it  is  perhaps  more  appropriate  to  specify 

stick  gain  in  units  of  volts/lb  rather  than  volts/inch.  We 

note  from  Table  5 that  feedthrough  accounted  for  an  error  score 
2 

of  0.009  volts  for  the  light  stiff  stick.  As  this  prediction  was 
obtained  for  a stick  gain  of  0,9  volts/lb,  we  may  now  restate 
equation  (18)  as 

a?  = .011  K?  (10) 
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where  K is  the  stick  electrical  gain  in  volts/lb.  This 

c 

relationship  is  shown  in  Figure  14  by  the  dotted  line. 

The  pilot/vehicle  model  is  used  to  predict  the  effect  of 
stick  gain  on  the  portion  of  the  tracking  error  variance  corre- 
lated with  the  tracking  input.  Except  for  the  stick  gain, 
problem  variables  and  pilot-related  parameters  are  as  described 
in  the  preceding  example  for  the  light-stiff  stick  configuration. 
The  relationship  between  input-related  error  score  and  stick  gain 
is  shown  by  the  dashed  curve  in  Figure  14. 

The  solid  curve  shown  in  Figure  14  shows  the  total  error 
variance  score  as  a function  of  stick  gain.  Minimum  error  variance 
is  predicted  for  a gain  of  about  6 volts/lb.  If  we  consider  as 
"optimal"  any  stick  gain  that  allows  an  error  score  within  10% 
of  the  theoretical  minimxim,  then  the  range  of  optimal  gain  extends 
between  approximately  0.2  and  1.3  volts/lb  for  the  particular 
tracking/biodynamic/vibration  situation  explored  in  this  example. 

The  range  of  optimal  stick  gains  predicted  here  is  consi- 
derably less  than  that  usually  associated  with  tracking  under 
static  conditions.  Although  stick  gain  was  not  explored  in  this 
study,  a recent  study  of  tracking  under  vibration  conditions 
confirms  the  prediction  that  considerably  less  lattitude  is 
allowed  in  selecting  stick  gains  for  vibration  environments  [8] . 

Theoretically,  analysis  via  the  pilot/vehicle  model  would 
be  sufficient  to  determine  the  range  of  acceptable  stick  gains  in 
the  absence  of  vibration  (or  in  vibration  situations  where  feed- 
through is  negligible).  In  this  situation,  however,  the  model 
would  have  to  be  augmented  to  account  for  the  upper  bound  on 
stick  gain  that  arises  from  neuromuscular  tremor  and  other  souces 


STICK  GAIN  (VOLTS/POUND) 


I'iquro  14.  rffect  of  Stick  Gain  on  Predicted  Trror  Variance 
(Light  Stiff  Stick) 


of  irreducible  random  control  activity.  We  suspect  that  a "residual 
motor  noise"  - akin  to  the  residual  observational  noise  parameters 
now  included  in  the  model  formulation  — would  serve  this  purpose. 
Additional  research  would  be  needed  to  quantify  this  model  para- 
meter. 

4.2.3  Effect  of  Stick  Damping 

In  this  final  example  we  use  the  design  procedure  to  explore 
the  effects  of  s.tick  damping  on  performance.  We  again  restrict 
analysis  to  the  light  stiff  stick  in  the  task  situation  defined 
in  the  initial  example,  and  the  results  are  presented  as  a design 
curve  showing  the  relation  between  error  variance  and  stick  damping. 

Analysis  proceeds  essentially  as  in  the  preceding  two 
examples.  Since  stick  damping  - and  thus  stick  impedance  — is 
a variable  of  the  analysis,  the  expression  of  equation 
predicts  the  relation  between  vibration-correlated  error  and 
stick  damping.  Except  for  the  variable  stick  damping,  all 
problem  parameters  are  the  same  as  defined  for  the  light-stiff 
configuration  in  the  initial  numerical  example.  The  effects  of 
stick  damping  on  vibration-correlated  error  score,  are  displayed 
by  the  dotted  curve  in  Figure  15. 

Model  analysis  allows  us  to  predict  the  effects  of  stick 
damping  on  the  tracking-related  component  of  the  error  variance 
score.  Except  for  stick  damping,  model  parameters  and  problem 
variables  are  as  defined  in  the  initial  example.  The  results  of 
this  analysis  are  displayed  by  the  dashed  curve. 
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The  effects  of  damping  on  totai  performance  is  shown  by 
the  solid  curve.  Because  of  the  small  contribution  of  feedthrough 
to  the  error  score,  only  a slight  reduction  in  score  can  be 
achieved  by  increasing  the  damping.  As  damping  is  increased 
beyond  about  2 lb/ (in/sec),  the  mechanical  impedance  of  the 
stick  begins  to  impair  tracking  efficiency,  causing  an  increase 
in  the  overall  error  variance  score. 


APPENDIX  A 

COKPUTATIONAL  TECHNIQUES 
A,1  Amplitude  Statistics 

standard  deviation  scores  were  computed  for  all  tracking 
and  biodynamic  variables.  These  scores  were  computed  over  a 
100-second  measurement  interval  that  commenced  approximately  15 
seconds  after  initiation  of  the  tracking  run.  As  all  variables 
had  essentially  zero  mean,  we  have  used  the  term  "rms  performance 
score"  elsewhere  in  this  report  to  signify  the  standard  deviation 
scores. 


Mean  performance  scores  were  computed  across  subjects,  and 
the  standard  deviation  (SD)  of  the  individual  scores  was  computed 
to  provide  a measure  of  subject- to-subject  variability.  Also 
computed  was  the  standard  error  (SE)  associated  with  the  mean 
to  provide  an  estimate  of  the  repeatibili ty  of  the  experiment. 

The  SD  and  SE  scores  were  computed  as  follows: 


SD 


1/2 


SE 


SD 

/N 


where  0^  is  the  average  performance  score  for  the  i^^  subject 
performing  a given  experimental  task,  o is  the  mean  of  the  o^, 
and  N is  the  number  of  test  subjects. 
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A.?  Power  Spectra 

Power  sp  ctra  were  computed  using  fast-Fourier  transform 
techniques  based  on  the  Cooley-Tukey  method  of  computing  trans- 
forms [12].  Spectra  were  com.puted  from  data  covering  a 100- 
second  measurement  interval  beginning  approximately  15  seconds 
after  the  initiation  of  the  tracking  run.  This  measurement 
interval  allowed  the  fundamental  analysis  frequency  of  about 
0.0628  rad/sec  to  be  an  integral  sub-multiple  of  the  theoretical 
base  frequency  (about  0.126  rad/sec'  used  in  constructing  the 
sums-of-sinusoids  tracking  input. 

The  procedure  for  differentiating  between  input-correlated 
and  remnant-related  components  of  a given  spectrum  was  similar 
to  that  used  in  the  past  ( 4 ) , except  that  correction  was  made 
for  non-ideal  inputs.  Let  us  first  describe  the  procedure  used 
with  ideal  inputs,  then  discuss  the  modifications  needed  in  this 
study. 


We  define  "ideal"  inputs  as  consisting  of  a number  of 
undistorted  sinusoids  that  are  harmonics  of  some  fundamental 
frequency.  Sinusoidal  components  of  the  tracking  input  and 
vibration  input  would  hiave  the  same  fundamental.  A measurement: 
interval  for  the  FFT  would  be  chosen  equal  to  the  reciprocal  of 
the  fundamental  input  frequency  (or  an  integral  multiple  thereof) . 
Under  these  circumstances,  signal  power  measured  at  a frequency 
not  coincident  with  a component  of  either  the  tracking  or  vibration 
input  would  be  defined  as  "remnant";  power  measured  at  a component 
frequency  of  the  tracking  input  would  be  assumed  to  be  correlated 
with  the  tracking  input  (provided  the  remnant  power  at  neighboring 
measurement  frequencies  was  sufficiently  small) ; and  power 
measured  at  a component  frequency  of  the  vibration  input  would  be 
considered  correlated  with  the  vibration  source  (again,  provided  a 
favorable  signal/noise  ratio) . 
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The  input  signals  obtained  from  the  data  tapes  were  non- 
ideal; that  is,  the  component  sinusoids  were  "smeared"  to  the 
extent  that  appreciable  signal  power  was  observed  over  two  or 
three  measurement  frequencies  on  either  side  of  the  nominal 
input  frequency.  Possible  sources  of  such  distortion  include 

(a)  imperfections  in  the  devices  used  to  generate  the  inputs  at 
experiment  time,  (b)  speed  fluctuations  of  the  analog  tape 
drive  during  recording  and/or  playback,  (c)  inaccuracies  or 
fluctuations  in  sampling  rate  during  digitization.  In  any 
event,  it  was  necessary  to  modify  the  analysis  program  so  that 
signal  power  appearing  at  measurement  frequencies  close  to  the 
nominal  input  frequency  would  be  properly  decomposed  into  input- 
correlated  and  remnant  power,  rather  than  being  considered 
entirely  as  remnant  power. 

The  correction  procedure  is  based  on  the  following  two 
assumptions;  (a)  signal  power  measured  at  a nominal  input 
frequency  is  not  significantly  corrupted  by  noise  and  can  be 
considered  to  be  entirely  correlated  with  the  input  signal; 

(b)  the  linear  transfer  characteristic  relating  the  response  of 
a particular  signal  to  the  input  is  essentially  invariant  over  a 
small  range  of  freq\iency  (say,  + 0.2  rad/sec).  As  we  show  below, 
the  first  assumption  can  be  tested  in  a given  measurement 
situation.  If  the  assumption  is  violated,  we  consider  the 
measurements  unreliable  and  proceed  accordingly.  We  cannot 
directly  test  the  second  assumption  in  a specific  situation. 
However,  analysis  of  previous  experimental  data  (including  model 
analysis)  indicates  that  describing  functions  and  spectra  vary 
smoothly  with  frequency  and  can  thus  be  considered  constant  over 
the  small  range  of  frequency  we  are  considering  here. 
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Let  Xj^  be  the  Fourier  transform  of  the  signal  x(t)  as 
computed  for  the  input  frequency.  According  to  the  first 
ass\uT\ption, 


^k  ^k  • ^k 

where  Ij^  is  the  transform  of  the  k^j^  input  component,  and  Tj^  is 
the  linear  transfer  characteristic  relating  the  two  signals.  At 
neighboring  frequencies,  where  both  input-correlated  and  remnant 
power  may  be  important,  the  following  relationship  applies: 


^k+j  ' ^ik+j 


(A-1) 


where  the  subscripts  "i"  and  "r”  denote  input-correlated  and  remnant- 
related  signal  components.  Since  the  transfer  characteristic  is 
assumed  invariant  for  small  values  of  the  index  "j",  the  correlated 
power  may  be  written  as 
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or,  in  terms  of  measurable  quantities; 
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Equations  (A-1)  and  (A-2)  may  be  combined  to  yield  the  following 
computation  for  the  remnant-related  component  of  the  signal: 
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We  now  define  the  "corrected"  measure  of  input-correlated  power  to 
be 
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and  the  remnant  power  associated  with  the  input  component  is 
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Since  pilot  remnant  appears  to  be  a continuous  function  of 
frequency,  the  estimate  of  remnant  power  is  normalized  by  the 
factor  2J  to  yield  average  power  per  measurement  interval.  No 
such  normalization  is  performed  on  the  estimate  of  input-correlated 
power,  as  the  theoretical  input  signal  has  power  at  only  a few 
widely-separated  frequencies. 

A. 3 Describing  Functions 

Estimates  of  pilot  describing  functions  at  input  frequencies 
were  obtained  by  dividing  the  transfom.  of  the  control  signal  by 
the  transform  of  the  error  signal.  Similar  manipulations  of 
transforms  were  performed  to  provide  estimates  (at  vibration 
frequencies)  of  the  transfer  characteristics  relating  (a)  shoulder 
acceleration  to  platform  acceleration,  (b)  inboard  bitebar  accelera 
tion  to  platform  acceleration,  (c)  outboard  bitebar  to  inboard 
bitebar  accelerations,  and  (d)  elbow  to  shoulder  accelerations. 


A. 4 Measurement  Reliability 

Measurement  reliability  is  determined  by  the  ratio  of  "signal^ 
power  to  "noise"  power  at  a given  measurement  frequency.  Two  signal/ 
noise  ratios  must  be  considered  when  analyzing  measures  of  tracking 
performance:  the  ratio  of  total  power  to  remnant  power  estimated  at 

a given  frequency,  and  the  ratio  of  remnant  power  to  system  (other 
than  pilot)  noise. 

As  noted  above,  reliable  estimates  of  the  input-correlated 
stick  and  error  power  are  needed  to  obtain  reliable  estimates  of 
the  pilot  describing  function.  The  best  estimate  of  the  input- 
correlated  power  at  a given  input  frequency  is  the  difference 
between  the  total  power  at  that  frequency  and  the  average  of  the 
remnant  pov/er  at  neighboring  frequencies.  There  is  some  un- 
certainty associated  with  the  estimate  of  input-correlated  power, 
however,  since,  at  a given  frequency,  the  remnant-related  signal 
may  combine  positively  or  negatively  with  the  input-correlated 
signal.  If  we  assume  that  the  phase  relation  between  remnant 
and  input-correlated  power  is  uniformly  distributed  between  0 
and  360  degrees,  the  standard  deviation  of  the  error  (in  dB) 
associated  with  the  estimate  of  input-correlated  power  is  related 
to  the  signal/noise  ratio  as  shown  in  Table  A-1.  In  this  case, 
"signal"  refers  to  input-correlated  power  and  "noise"  to  remnant 
power. 


Since  vibration  influences  tracking  performance  primarily 
through  an  increase  in  pilot  remnant,  reliable  estimates  of  the 
remnant  spectrum  are  desired  to  test  models  of  pilot  performance 
in  vibration  environments.  Table  A-1  can  be  used  to  determine 
the  reliability  of  the  remnant  estimate  as  well,  where  we  now 
interpret  the  remnant  as  "signal"  and  system  noise  as  "noise". 


Note  that  system  noise  can  only  be  obtained  in  an  experiment  in 
which  there  is  no  tracking,  or  in  which  tracking  is  performed  by 
an  appropriate  electronic  feedback  network. 

For  the  data  analyzed  in  this  study  we  adopted  a signal/ 
noise  criterion  of  3 dB,  which  corresponds  to  approximately  a 
7 dB  ratio  of  total  power  to  noise  power.  That  is,  spectral  and 
describing  function  measurements  were  considered  unreliable  and 
disregarded  in  further  analysis  wherever  total  signal  power  failed 
to  be  at  least  7 dB  greater  than  the  estimated  noise  pov;cr. 
"Corrected"  (as  described  above)  estimates  of  remnant  power  were 
used  in  making  this  test. 

Vibration  spectra  from  one  experimental  run  were  examined 
in  considerable  detail  to  determine  the  fidelity  of  the  platform 
vibration  spectrum  and  to  assess  the  utility  of  deriving  linear 
models  for  biodynamic  response  mechanisms.  The  spectrum  of  the 
platform  acceleration  was  examined  in  the  vicinity  of  the 
nominal  vibration  frequencies  to  explore  for  possible  distortion 
in  the  component  sinusoids.  The  spectrum  appeared  to  be  very 
clean  for  tlie  particular  data  examined.  Power  was  down  by  at 
least  20  dli  outside  a band  of  +0.2  rad/sec  about  cacii  nominal 
input  frequency.  Furthermore,  there  were  no  significantly  large 
"spiJees"  of  power  Ijrtween  input  frequencies.  Thus,  power  at  non- 
input frequencies  was  relative) y smoothly  distributed  over 
frequency,  at  least  in  the  frequency  reoion  spanned  Jjy  the  five 
viln-ation  input  f roquenci<.s . 
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Detailed  Inspection  of  the  spectra  for  shoulder  and 
bitebar  signals  showed  that  the  fall-off  about  the  nominal  vibration 
frequency  was  about  the  same  as  the  fall-off  of  the  platform 
spectrum.  Furthermore,  the  ratio  of  vibration-correlated  to 
remnant  power  (measured  in  the  vicinity  of  vibration  input 
frequencies)  was  generally  quite  large  for  all  biodynamic  quan- 
tities. 


Biodynamic  data  from  selected  experiments  were 
analyzed  to  determine  the  fraction  of  total  response  power 
not  correlated  with  the  platform  acceleration.  As  shown  in 
Table  A-2,  the  fractional  remnant  power  was,  on  the  average, 
about  20%.  Head  and  limb  motions  were  thus  largely  due  to 
linear  coupling  with  the  platform,  and  emphasis  was  accordingly 
placed  on  measuring  biodynamic  describing  functions. 

A. 5 Equivalent  Continuous  Spectra 

Tl\e  100-second  measurement  interval  yielded  line  spectra 
with  a spacing  of  approximately  0.0C23  rad/nec  between  successive 
lines.  Those  line  spectra  were  transf omied  into  equivalent 
continuous  spectra  Jor  comyjarison  with  model  results.  Remnant- 
related  spectra  were  transformed  simply  by  dividing  the  power  at 
ouch  measurement  frequency  l)y  the  ”measurem€‘nt  window"  of  0,0623 
rad/soc . 


Because  the  power  in  the  tracking  input  was  basically 
concentrated  at  five  frequencies  instead  of  being  distributed 
continuous] y over  irequoncy,  input-correlated  spectral  measure- 
ments could  not  be  transformed  into  equivalent  continuous  density 
functions  by  normalizing  witli  respect  to  the  measurement  sindow. 
Instead,  a "box  car"  approximation  was  made  to  the  input  spectrum. 
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That  is,  each  sinusoidal  component  of  the  input  was  considered  to 
represent  the  integral  of  a constant  power  density  extending  over 
the  region  of  fr  ,quency  space  bounded  by  the  geometric  midpoints 
of  the  component  frequency  and  the  two  neighboring  frequencies. 
Input*  correlated  power  measurements  were  then  normalized  with 
respect  to  the  frequency  interval  spanned  by  the  corresponding 
segment  of  the  box-car  input. 

A. 6 Some  Definitions 

We  use  the  term  "spectrum"  to  refer  both  to  discrete 
spectral  measurements  as  well  as  to  continuous  power  spectral 
density  functions.  Spectra  are  ciefined  as  beiiig  nonzero  over 
positive  frequencies  only,  and  signal  variance  is  related  to 
the  "spectruni"  as  follows: 

OO 

2 

0 

2 

where  is  the  spectrum  of  the  signal  x(t)  and  is  the 

variance  of  x(t).  In  the  case  of  a discrete  spectrum,  the 
integration  over  fr'quency  is  to  be  interpreted  as  a summation. 

Tracking  inputs  and  vibration  inputs  are  assumed  to  be 
linearly  uncorrelated.  By  definition,  "remnant"  power  is  not 
linearly  correlated  with  either  of  these  signals.  Therefore, 
the  spectrum  of  any  tracking  or  biodynamic  variable  may  be 
considered  to  be  the  sum  of  three  component  spectra: 

■ (t)  (1)  4*  ^ 

XX  X^X^  -f  X^X^  X^Xj, 
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where  the  subscripts  "i",  ''v",  and  "r",  designate  components  of 
signal  power  that  are  correlated  with  the  tracking  input, 
correlated  with  the  vibration  input,  and  remnant-related. 
Similarly,  signal  variance  may  be  considered  as  the  combination 
of  three  components: 
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where  each  component  variance  is  equal  to  the  integral  over 
positive  frequencies  of  the  corresponding  component  spectrum. 
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APPENDIX  B 

ANALYSIS  OF  BIODYNAMIC  RESPONSE  DATA 

There  are  a number  of  blodynamic  response  mechanisms  that 
act  to  degrade  tracking  performance  when  the  pilot/vehicle  system 
is  subjected  to  vibration  environments.  Direct  mechanical 
coupling  between  the  vibration  input  and  the  control  stick 
(i.e.i  vibration  feedthrough)  will  introduce  a control  input 
that  is  linearly  correlated  with  the  system.  Such  an  input  can 
have  a significant  effect  on  performance,  especially  if  the 
vehicle  has  resonant  characteristics  at  vibration  frequencies 
(e.g.,  aircraft  bending  modes).  Feedthrough  may  also  increase 
the  level  of  wide-bond  pilot  remnant  by  interfering  with  neuro- 
muscular feedbacks  needed  for  precise  control  of  the  stick. 
Similarly,  vibration  of  the  limbs  may  also  interfere  with  these 
feedbacks . 

Relative  motion  between  the  eye  point-of-regard  and  the 
display  may  induce  blurring  and  thus  degrade  visual  inputs  needed 
for  effective  control.  Since  the  oculomotor  sys<  "^m  will  not  allow 
the  eye  to  compensate  for  relative  motions  at  frequencies  above 
4 Hz,  visual  interference  effects  can  bo  inferred  to  a large 
e>:tent  from  a study  of  vibration-induced  head  motion.  Finally, 
exposurj  of  long  duration  to  substantial  vibration  can  be 
i.,xp(>cted  to  induce  fatigue  and  and  a general  feeling  of  ill-being, 
thereby  leading  to  further  degradation  in  tracking  performance. 

A comprehensive  set  of  biodynamic  response  measures  was 
obtained  in  the  experimental  program  so  that  the  relative  importance 
of  various  vibration-related  interference  mechanisms  could  bo 
assessed.  Acceleration  of  the  vibrating  platform  was  recorded 
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along  with  shoulder  and  elbow  acceleration.  Accelerometers  mounted 
on  a bitebar  provided  measures  indicative  of  linear  and  rotational 
head  motion.  In  addition#  the  Fourier  analysis  techniques  described 
in  Appendix  A allowed  us  to  identify  the  portion  of  the  pilot's 
control  response  linearly  related  to  the  platform  vibration.  From 
these  measures  were  obtained  linear  models  (i.e.#  describing 
functions)  relating  body#  head#  and  limb  motions  to  the  vibration 
input.  These  models  were  then  used  to  infer  measures  that  were  not 
directly  obtained  (such  as  relative  motion  between  the  eye  fixation 
point  and  the  display). 

Response  measures  obtained  for  the  various  biodynamic 
mechanisms  are  presented  below#  and  special  analysis  procedures 
used  in  obtaining  these  meeisures  are  described.  We  begin  presen- 
tation of  the  experimental  results  with  a description  of  the  data 
base , 


B , 1 , 1 Data  Base  for  Biodynamic  Measurements 

As  noted  in  Section  2#  not  all  test  subjects  participating 
in  the  experimental  program  participated  in  each  sub-experiment. 

In  addition#  as  explained  later  in  this  Appendix#  certain  bio- 
dynamic measures  obtained  in  the  early  experiments  could  not  be 
calibrated  reMably  for  analysis. 

The  data  base  for  the  results  presented  in  this  Appendix  is 
indicated  in  Table  D-1,  An  "X”  in  the  column  labeled  C (or  S) 
indicates  that  a reliable  measure  was  obtained  with  the  stic)( 
located  in  the  center  (or  side)  position. 


Light  Spring 

X 
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Heavy  Spring 

Light  Stiff 
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B , 1 . 2 Rms  Acceleration  Scores 

As  noted  in  Tabic  B-1/  only  the  light  spring  and  light  stiff 
stick  configurations  yielded  reliable  measures  of  all  biodynamic 
varia’  les  for  both  center  and  side  locations,  Rms  acceleration 
scr  as  for  shoulder,  elbow,  inboard  bitebar,  and  outboard  bitebar 
accelerometers  are  shown  for  both  locations  in  Table  B>2a,  Scores 
(in  units  of  g)  are  shown  for  each  of  five  test  subjects,  and  the 
mean  performance  of  these  subjects  is  given. 

A summary  of  the  analysis  of  variance  of  these  scores  is 
given  in  Table  B-2b.  Entries  show  the  level  of  significance  of 
performance  differences.  No  significant  primary  effects  due  to 
st.’.ck  location  or  stick  configuration  Are  found  for  any  of  the 
four  biodynamic  measures,  and  significant  differences  across 
subjects  are  shown  only  for  the  bitebar  accelerations.  Since 
stick  location  has  no  significant  effect  on  these  measures, 
further  analysis  of  biodynamic  data  is  based  ( /here  possible) 
on  measures  averaged  across  side  and  center  stick  locations. 

Rms  acceleration  scores  are  shown  for  seven  test  subjects 
in  Table  B-3.  Also  shown  are  mean  performance  as  well  as  standard 
deviations  and  standard  errors  computed  acorss  subjects.’* 

Wherever  possible,  data  have  been  averaged  across  side  and  center 
stick  locations. 


*^See  Appendix  A for  a definition  of  statistical  measures. 


Table  B~2 

Effect  of  Stick  Location  on  Rms  Acceleration  Scores 


PH  DM 


a.  Shoulder 

iLighv.  Spring  Center 

Light  Spring  Side  .23 

Light  Stiff  Center  .30 

Light  Stiff  Side  .240 


307  .275 

254  .308 


.319 

.267 


. E Ibow 


Light  Spring 
Light  Spring 
Light  Stiff 
Light  Stiff 


.282 
. 36  3 
.708 


7 .236 

4 I .327 
6 ' .380 


.402 

.338 

.421 


Inboard  Ditebar 


Center 


367 

. 365 

391 

. 36  2 

346 

. 347 

Light  Spring 


Center 
de 
nt 
de 


.311 
. 325 
. 306 


.3 

.3 
.274 


Variable 

Rlbow 

Inboard  Ditebar 
Outboard  Ditebar 


ource  of  Difference 
1 [tocation  ILocatioh 

Location  Stick  i Subject  X X 

Stick  Subject 


X X 

Subject  Subject 


Entries  show  level  of  significance 
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Table  B-3 


Rms  Acceleration  Scores 


’ Stick 

"1 

Individual  Subjects 

Statistics 

>—  i>.  II  ■^l^  ■ 1. 

CBH 

PH 

RM  i DM  JP  1 AP 

d“^ 

Mean 

SD  SE 

a.  Shoulder 


; Light  Spring  i 

T?3T 

' T2  98‘ 

TT751 

mrr 

TTeTTTOT* 

THT 

T5?T 

“Tm‘r&Tg 

1 Medium  Spring: 

. 343 

.275 

.323 

.313j .290 

. 233 

.296 

.039  .016 

Heavy  Spring 

. 276 

. 260 

.272 

.194 

.256; .298 

.269 

.261 

.032  „012 

, Light  Stiff  i 

. 270 

.280 

.200 

.292 

.394 : .276 

. 298 

.287 

.057  .022 

Medium  Stiff 

. 286 

.314 

.301 

.326 

.362; .268 

.239 

.299 

.040  .015 

Heavy  Stiff  i 

.288 

.317 

.257 

. 286 

. 233 i .286 

. 282 

.278 

.027  .010 

b.  Elbov/ 


Light  Spring  \ 

.605^ 

.306 

.269' .296 

.308  .287 

. 278 

. 335 

.119. .045 

Medium  Spring! 

. 794  , 

. 780 

' .414 

.326  .220 

. 428 

.494 

.239 

.098 

Heavy  Spring  j 

. 302  : 

. 353 

.323' . 837 

.365; .524 

. 347 

.436 

.191 

.072 

Light  Stiff  ! 

. 342  ■ 

. 536 

.269-'  . 37  5 

.318  .354 

.290 

. 355 

.088 

.033 

Medium  'itiff  ! 

. 354 

. 353 

.6961 . 408 

. 336( . 256 

. 424 

.404 

.140 

. 053 

Heavy  Stiff 

.415 

. 435 

. 364! . 374 ! . 388!  .366 

. 376 

. 388 

.027 

c.  Inboard  Ili  tebar 


Light  Spring 

. 326 

. 252 

. 346 

. 354  . 280 

. 380 

. 31  5 

. 322; 

.044 

. 017 

Medium  Spring 

. 33  5 

. 292 

.360  .275 

. 352 

.422 

. 339' 

.053 

. 022 

Heavy  Spring 

0 342 

. 296 

. 372 

.327  .302 

. 324 

. 29  5 

. 32  3 

. 02  0 

. Oil 

Light  Stiff 

. 397 

.247 

. 283 

.368  .255 

. 354 

.312 

.317 : 

.058 

.022 

Medium  Stiff 

. 322 

. 274 

.450 

.388  .280 

.361 

. 369 

. 349 

. 06  2 

. 024 

Heavy  Stiff 

. 354 

.275 

.344 

. 308, , 345 

. 369 

. 27,3 

.32  4 

. 03  9 

. 0 :i,  5 

! 

d.  Outboard 

Hi  teliar 

Light  Spring 

. 288- 

. 220 

. 300 

.2'9?  75irr 

. W 

. 28  IT.  2 83 

• 

. 012 

Medium  Spring 

. 312 

.277 

.317  .256 

. 315 

. 368 

. 30  8-. 

.038 

. 016 

Heavy  Spring 

.165 

. 227 

. 246 

.212  .212 

.406 

. 30  1 

. 25  3' 

. 079 

. 0 30 

Light  Stiff 

. 316 

. 232 

.238 

.305  .242 

. 308 

. 27  8 

.274 

. 036 

. 014 

Medium  Stiff 

. 306 

.267 

. 494 

.330  .261 

.323 

. 327 

. 330 

. 07  8 

. 029 

Heavy  Stiff 

. 195 

.285 

. 249 

. 245.  . 237 

.492 

. 195 

. 271. 

.102 

.0  39 

All  scores  are  in  q's. 
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Table  B-4  contains  a summary  of  the  analysis  of  variance 
performed  on  the  scores  contained  in  Table  B-3.  Because  of  the 
missing  data  point  for  Subject  RM,  this  analysis  is  performed 
on  data  provided  by  the  remaining  six  subjects.  Since  there  is 
only  one  entry  per  subject  per  stick  configuration  (center,  side, 
or  the  average  of  the  two),  we  can  compute  only  main  effects. 

Table  B-4  shows  that  there  are  no  main  effects  for  the 
shoulder  and  elbow  accelerations  and  that  only  inter-subject 
differences  are  significant  for  the  bitebar  accelerations. 

We  therefore  conclude  that  limb  and  head  vibrations  are  essen- 
tially invariant  with  regard  to  the  specifics  of  the  control 
stick.  Furthermore,  the  lack  of  inter-subject  differences 
(at  least  as  far  as  shoulder  and  elbow  response  is  concerned) 
justifies  consideration  of  the  "average  man"  that  is  defined 
by  the  average  response  behavior  of  the  test  subjects, 

B.2  Vibration  Feedthrough 

B , 2 . 1 Analysis  Procedures 

Vibration  feodthrough  is  defined  as  the  component  of  the 
pilot's  control  response  that  is  linearly  related  to  the  vibration 
input.  Describing  functions  relating  control  response  to  vibra- 
tion were  obtained  using  the  Fourier  analysis  procedure  described 
earlier.  These  measures  were  averaged  across  subjects  for  pur- 
poses of  model  analysis.  In  order  that  the  effects  of  stick 
location  be  observable,  however,  av’eraging  was  n<^  performed 
across  center  and  side  conditions. 
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Table  B-4 

Analysis  of  Variance  of  Data  in  Table  B-3 


Source  of  Difference 


Variable 


Shoulder 

1 

Elbow  i 

Inboard  Bitebar 

.001 

Outboard  Bitebar 

1 

j 

i .01 

1 

Entries  show  level  of  significance 
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As  expected f the  feedthrough  describing  function  varied 
with  the  control-stick  characteristics.  In  order  to  characterize 
these  results  in  a manner  that  would  allow  extrapolation  to 
control-stick  designs  not  studied  in  this  program,  a regression 
analysis  was  performed  to  identify  portions  of  the  response 
mechanism  independent  of  control-stick  parameters. 


The  following  model  structure  was  assumed: 


where 


- rgTflr^kTTi.T  <®-i> 

C « Stick  response  in  inches 
a a table  vibration  in  g's 

ZT  ■ the  effective  biodyramic  transfer  impedance  of  the 
pilot  in  Ib/g, 

ZO  = the  effective  biodynamic  output  impedance  of  the 
pilot  in  Ib/in, 

ZS  = the  mechanical  impedance  of  the  control  stick  in 
Ib/in. 


The  parameters  ZT  and  ZO  are  intended  to  reflect  the  biodynamical 
properties  of  the  pilot  and  the  interface  between  the  pilot  and  the 
cockpit.  For  this  model  form  to  have  maximum  predictive  capability, 
ZT  and  ZO  should  be  independent  of  the  mechanical  and  electrical 
properties  of  the- control  stick.  (Location  of  the  control  stick 
will  in  general  affect  these  parameters,  however,  since  stick 
location  influences  pobitioning  of  the  body  and  limbs.) 


Two  sets  of  transfer  characteristics  were  determined:  one 

set  for  center  stick  location  and  one  for  side.  A non-parametric 
regression  analysis  was  performed,  that  is,  no  analytic  form  was 
assumed  for  the  frequency  dependency  of  ZT  and  ZO.  Instead, 
independent  regression  analysis  Wcib  performed  on  the  data  obtained 
at  each  vibration  frequency. 
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The  algorithm  for  the  regression  analysis  is  developed  as 
follows.  Let  Xj^(jo)j^)  be  the  stick/table  transfer  function 
measured  at  a particular  frequency  with  the  i^^  control 
configuration,  and  let  X^(jwj^)  be  the  transfer  that  is  pre- 
dicted by  the  model  of  equation  (B-1)  for  the  "best"  values  of 
ZT  and  ZO.  From  equation  (n-1)  we  obtain 


A 


ZT 

ZO  + ZS^ 


(B-2a) 


i-  = i^rtZO  + ZS^]  (B-2b) 

^i 


where  the  argument  (jw^)  is  omitted  for  compactness. 


We  now  define  ZTI  = 1/ZT 
ZOT  = ZO/ZT 


(B-3) 


which  allows  equation  (B-2b)  to  be  rewritten  as 

= ZOT  + ZTI  • ZS. 

^ 1 

^i 

The  fractional  error  in  1/X.  (which  is  approximately  equal  in 

^ A 

magnitude  to  the  fractional  error  in  for  small  errors)  is 


1 1 


= ZOT  • X^  + ZTI  • Z^  . ZS^  - 1 
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The  expression  for  the  modeling  error  may  be  given  in  the 
following  matrix  form 


which  we  may  express  by  the  shorthand  notation 


E ■ A • 2 - B 


Note  that  the  matrix  A consists  of  measured  (X^)  and  known  (ZS^^) 
quantities,  both  of  which  will  generally  be  complex. 


ZOT  and  ZTI  can  je  determined  if  measurements  are  obtained 
for  two  or  more  stick  configurations.  The  values  that  minimize 
the  sum  of  the  squared  modeling  errors  are  given  by  the  expression 


(A*  A) 


* 

A'  B 


where  A'  denotes  the  complex  conjugate  of  the  transposed  matrix. 
The  parameters  ZT  and  ZO  are  then  derived  from  the  identities  of 
equation  (B-3). 

If  measurements  are  obtained  for  only  two  stick  configura- 
tions, values  for  ZT  and  ZO  can  be  found  tr  provide  a perfect 
match  to  the  stick/table  transfer  function  at  each  measurement 
frequency.  If  measurements  are  obtained  for  more  than  two  stick 
configurations,  a good  match  to  the  data  will  be  possible  only  if 
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the  basic  assumption  underlying  the  regression  analysis  is  valid; 
namely,  that  there  are  identifiable  components  of  the  vibration 
feedthrough  mechanism  that  are  independent  of  control-stick  para- 
meters. 


B.1.2  Experimental  Results 


Average  stick/platform  transfer  functions  are  given  in 
Figures  B-1  and  B-2  for  each  of  the  six  control  stick  configurations 
Measures  for  both  center  and  side  locations,  where  available,  are 
shown  separately*  Amplitude  ratios  (expressed  in  dE)  represr nt  the 
vibration-correlated  control  response  in  inches  per  g of  platform 
acceleration. 


The  amplitude  ’'atios  reveal  no  consistent  differences 
between  center  and  side  locations.  Location  does  have  a noticeable 
effect  on  phase  shift,  however.  A less  negative  phase  shift  is 
consistently  shown  for  the  center  location,  especially  at  the 
lower  vibration  frec;xu.-ncies.  In  addition,  the  curves  for  the 
side  location  show  an  increase  in  phase  shift  occurring  between 
the  5 and  7 Hz  (31  and  44  rad/sec)  vibration  ..  iguencies  — a 
trend  not  revealed  for  the  center  location.* 


The  transfer  and  output  impedances  derived  from  these 
measurements  are  shown  in  Figure  13-3.  Also  shown  for  comparison 
are  the  mechanical  impedances  for  each  of  the  six  control  sticks. 
Amplitude  ratios  indicate  pounds/g  for  the  transfer  impedance 
and  pounds/inch  for  the  output  and  stick  impedances.  Again, 


fsThce  phase  shift  is  a circular  measure,  it  is  possible  that  some 
of  the  phase  measures  shown  in  these  figures  are  in  error  by  360* , 
Additional  measures  in  the  region  of  5-7  Hz  (or  alternatively,  a 
relieible  theoretical  model)  are  required  to  shed  further  light  on 
the  frequency  dependency  of  the  phase  shift  response. 
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stick  location  has  an  appreciable  effect  on  phase  shift  and  a 
negligible  effect  on  amplitude  ratio.. 

The  impedance  functions  help  explain  the  way  in  which  the 
stick/platform  transfers  (particularly  the  amplitude  ratio)  change 
with  stick  parameters.  For  stick  configurations  which  yield  a 
stick  impedance  less  than  the  output  impedance  ZO,  one  would 
expect  the  stick/platform  transfer  to  be  largely  independent  of 
stick  impedance.  Conversely#  as  the  stick  impedance  become  much 
larger  than  ZO#  the  stick/platform  transfer  should  begin  to  vary 
inversely  with  stick  impedance. 

The  transfers  shown  in  Figures  B-1  and  B-2  largely  confirm 
these  predicted  trends.  The  light  spring  and  medium  spring  sticks 
show  similar  stick /platform  transfers.  On  the  other  hand#  the 
stiff  sticks  yield  stick/platform  amplitude  ratios  that  appear  to 
decrease  inversely  with  the  stick  impedance. 


Considerable  effort  was  made  to  find  simple  analytic 
approximations  to  the  transfer  and  output  impedances.  While  good 
second-order  fits  were  obtained  for  the  amplitude  ratios  for  both 
side  and  center  locations#  a good  match  to  the  phase-shift  behavior 
was  not  obtained.  Accordingly#  we  have  refrained  from  presenting 
smooth  approximations  to  these  impedance  functions. 


In  order  to  assess  the  validity  of  summarizing  vibration 
feedthrough  mechanisms  in  terms  of  the  impedance  functions  ZT  and 
ZO#  we  have  used  these  derived  measures  to  "predict"  the  vibration- 
correlated  component  of  the  control  scores  by  carrying  out  the 
following  computation: 


(volts 


5 

I 

i«l 


ZT 
ZO  > 


e 


Opttp 
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These  predicted  scores  are  compared  with  the  corresponding  experi- 
mental measurements  in  Figure  B-4.  For  the  most  part,  there  is  good 
agreement  between  predicted  and  measured  vibration-correlated  scores 
over  a large  dynamic  range.  Thus,  we  are  justified  in  using  the 
impedance  concept;  we  can,  with  reasonable  confidence,  extrapolate 
the  results  of  this  experiment  to  situations  involving  other  control 
configurations,  provided  we  restrict  the  problem  situation  to  one 
having  the  same  biodynamic  and  vibration  environments  explored  in 
this  study. 

B.3  Shoulder  and  Elbow  Response 
B.3.1  Analysis  Procedures 

Describing  functions  relating  shoulder  and  elbow  acceleration 
to  platform  accelerat.'.on  were  computed  as  described  in  Appendix  A. 
Because  of  the  constraints  of  the  computer  program  that  was  used  to 
digitize  the  analog  data,  however,  special  effort  was  required  to 
calibrate  acceleration  variables.  Since  reliable  calibrations 
could  not  be  obtained  for  all  experiments,  some  biodynamic  data  had 
to  be  omitted  from  analysis  because  of  scale-factor  uncertainties* 
Calibration  methods  are  discussed  below. 

The  computer  program  that  was  available  for  digitizing  the 
experimental  data  required  that  two  different  but  constant  signal 
levels  be  provided  for  each  analog  variable.  Given  these  two 
signal  levels,  and  provided  with  the  interpretation  of  each  signal 
level  in  terms  of  experimental  units,  the  program  was  then  able  to 
compute  for  each  problem  variable  the  linear  relationship  between 
signal  level  in  experimental  units  and  volts  received  at  the 
A/D  converter. 
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Figure  D-4,  Compari.son  of  Predicted  and  Measured  Vibration- 
Correlated  Control  Scores 
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Constant  levels  were  recorded  on  analog  tape  for  each  of 
the  tracking  variables,  and  calibration  of  these  signals  proceeded 
without  difficulty.  Since  the  accelerometers  could  not  be  subjected 
to  a constant  acceleration  for  a sufficient  length  of  time,  however, 
it  was  not  possible  to  record  the  desired  constant  cal..bration 
voltages  on  analog  tape.  Instead,  all  accelerometers  were  sub- 
jected to  sinusoidal  vibration  at  10  Hz,  and  this  information  was 
used  during  playback  to  generate  constant  calibration  levels. 

Calibration  of  acceleration  signals  was  considered  unreliable 
for  Experiments  1 and  2 because  (a)  calibration  signals  were  not 
"clean",  and  (b)  transformation  of  the  sinusoidal  calibration 
signal  to  a d.c.  level  was  not  sufficiently  sophisticated.  For 
these  experiments,  we  simply  examined  the  sinusoidal  calibration 
signal  visually  by  playback  into  an  oscilloscope,  and  used  the 
analog  computer  to  generate  d.c.  voltages  equal  to  plus  and  minus 
the  zero-to-peak  amplitude  of  the  sinusoid.  Accuracy  of  this 
procedure  was  hindered  by  the  presence  of  significant  low-frequency 
noise  on  the  10  Hz  sinusoid  that  was  recorded  on  analog  tape,  and 
the  short  duration  of  the  calibration  signal  prohibited  much  in 
the  way  of  electronic  filtering.  We  therefore  felt  that  accelera- 
tion signals  obtained  in  these  two  experiments  were  subject  to 
possible  significant  errors  with  regard  to  determining  the  scale 
factor  and  should  not  be  used  for  analysis. 

The  reliability  of  calibration  was  greatly  improved  in 
subsequent  experiments  by  providing  signals  less  corrupted  by 
noise  and  of  longer  duration,  and  by  using  an  electronic 
filtering  technique  to  determine  the  required  constant  calibra- 
tion levels.  Each  calibration  sinusoid  was  processed  by  a low- 
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pass  filter  to  reduce  the  possible  low-frequency  noise  corruption 
(including  d-c  bias) , and  the  filter  output  was  rectified  to 
yield  a voltage  equal  to  the  zero-peak  swing  of  the  sinusoid. 

B.3,2  Experimental  Results 

Describing  functions  relating  shoulder  acceleration  to 
platform  acceleration  are  shown  for  the  six  stick  configurations 
in  Figure  B-5.  These  results  reflect  measures  obtained  from  the 
center,  side,  or  average  of  center  and  side  locations,  according 
to  the  availability  of  the  data  (see  Table  B-1) . Since  no 
consistent  differences  are  observed  across  stick  configuration, 
we  are  justified  in  coalescing  these  data  into  a single  average 
describing  function,  also  shown  in  the  same  figure.  The  smooth 
curve  represents  an  analytic  approximation  to  the  average  shoulder/ 
platform  describing  function  having  the  form 

,!s(s)  = 1 + (B/K)s 

“‘p  1 + (B/K)s  + (|)s^ 

2 

Where  B/K  = .033  seconds  and  M/K  = .0012  seconds  . 

Elbow/shoulder  describing  functions  are  shown  for  all  stick 
configurations  in  Figure  B-6.  Except  perhaps  for  some  high- 
frequency  phase  differences  associated  with  t^e  stiff  sticks,  the 
elbow  appears  co  respond  in  phase  with  the  shoulder  but  at  a 
slightly  greater  amplitude.  By  averaging  across  the  five 
frequency  points  and  across  the  six  sticks  we  find  that  the 
amplitude  of  the  elbow  response  is  about  1.4  dB  greater  than  the 
shoulder  response.  This  is  a somewhat  lesser  difference  than 
indicated  by  the  rms  acceleration  scores  shown  in  Table  B-3, 
where  the  average  ratio  of  rms  elbow  score  to  rms  shoulder  score 
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i is  about  1.4  (a  difference  of  about  3 dB),  The  data  of  Figure 

) 

I B-6  are  based  entirely  on  the  components  of  shoulder  and  elbow 
I response  that  are  linearly  correlated  with  the  platform  vibration, 
j As  shown  previously  in  Table  A-2,  a larger  fraction  of  the  elbow 
I response  is  uncorrelated  with  platform  motion,  which  may  account 

I for  the  greater  difference  in  rms  acceleration  scores. 

) 

t 

i 

I 6.4  Head  Motion 

I Measurements  obtained  from  the  inboard  and  outboard  bitebar 

I accelerometers  were  used  to  compute  describing  functions  relating 
I translational  and  rotational  response  of  the  head  to  platform 
j acceleration.  Relative  motion  between  the  head  point“Of~regard 
and  the  display  was  also  computed  from  these  measures.  Computa- 
tional techniques  and  experimental  results  are  described  below. 

t 

I 

Analysis  Procedures 

The  measurement  situation  is  diagrammed  in  Figure  B-7. 

Let  D = the  viewing  distance,  inches 

S = separation  between  the  inboard  and  outboard  accelero- 
meters, inches 

- vertical  displacement  of  the  inboard  bitebar,  inches 

= vertical  displacement  of  the  outboard  bitebar,  inches 

= vertical  displacement  of  the  head  point-of-regard , inches, 
obtained  by  projecting  the  line  of  the  two  accelerometers 
to  the  plane  of  the  display 

9 = angle  of  head  rotation,  radians 
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For  small  angles  of  rotation#  the  following  approximations  are 
valid; 


(B-4) 


Oq  « 386.4 


<“o  ■ “1> 


. 386,4 


<“h  ~ °l’ 
D 


(B-5) 


where  Oj  # a , and  <Xu  denote  vertical  accelerations  of  the  inboard 
i'  o'  h 

accelerometer,  outboard  accelerometer,  and  the  head  point-of-regard 

2 

in  g;  and  Qq  denotes  rotational  acceleration  of  the  head  in  rad/sec 


Measurements  at  our  disposal  were  the  outputs  of  the  inboard 
and  outboard  accelerometers.  These  signals  were  Fourier  transformed 

ai 

to  provide  the  following  describing  functions:  (a)  relating 

acceleration  of  the  inboard  accelerometer  to  platform  accelerometer, 
and  (b)  relating  acceleration  of  the  outboard  accelerometer  to 
that  of  the  inboard  accelerometer.  Both  these  quantities  are 
dimensionless. 


Translation  of  the  head  was  considered  equivalent  to  trans- 
lation of  the  inboard  bitebar,  and  the  relationship  between  head 
rotation  and  platform  acceleration  was  computed  in  terms  of 
measureable  quantities  by  appropriate  manipulation  of  equation 
(B-Ei).  Thus, 


(B-6) 


386.4 


(B-7) 
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^0  2 

where  — is  a dimensionless  quantity  and  has  units  of  (rad/sec  )/g. 
P P 

Relative  motion  between  the  head  point-of-regard  and  the 
display  was  computed  as  follows.  First,  we  manipulate  the  identity 
shown  in  equation  (B-5)  to  yield  the  following  expression  for  head 
point-of-regard  acoeleration: 


a 


h 


) 


We  then  compute  the  describing  function  relating  acceleration  of 
the  head  point-of-regard  to  platform  acceleration  in  terms  of 
measureable  quantities; 


Finally,  we  note  that  the  display  is  attached  rigidly  to  the  plat 
form,  and  the  spectrum  acceleration  of  the  head  point-of-regard 
with  respect  to  the  display  is  given  as 


“dh^^dh 


Ui 


cx. 


- 1 


P P 
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^ - 1 

1^-1 
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i“i 

J“p 

P P 


(B-8) 


where  is  the  relative  acceleration  of  the  head  point-of- 

regard  in  inches/sec^  and  'J’apUp  is  the  platform  acceleration  in 
g.  Appropriate  integrations  of  this  spectral  quantity  yield 
mean-squared  displacement  and  velocity  of  the  head  point-of-regard 
relative  to  the  display. 
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B.4.2  Experimental  Results 


Describing  functions  relating  inboard  bitebar  acceleration 
to  platform  acceleration  are  shown  in  Figure  B-8  for  the  six  stick 
configurations;  describing  functions  relating  outboard  to  inboard 
bitebar  accelerations  are  shown  in  Figure  B-9*  Results  obtained 
from  side  and  center  locations  have  been  combined  where  possible. 

The  bitebar/platform  describing  functions  are  virtually 
identical  for  the  six  stick  configurations;  hence,  we  may  consider 
the  average  describing  function  (also  shown  in  Figure  B-8) . 

Figure  B-9  shows  that  the  describing  functions  relating  outboard 
to  inboard  bitebar  acceleration  are  nearly  identical  for  four 
of  the  stick  configurations  (all  but  heavy  spring  and  heavy 
stiff).  These  transfer  functions  are  essentially  flat  with 
near-unity  gain  and  are  averaged  together  in  Figure  B-9c.  The 
remaining  describing  functions  show  an  amplitude  variation  of 
about  4 dB  across  the  measurement  range.  It  is  not  clear  why 
these  transfers  should  differ  from  the  others,  and  in  Section 
B.4.3  we  discuss  the  possibility  of  measurement  error. 

Describing  functions  relating  head  rotational  acceleration 
to  platform  acceleration  were  computed  from  the  preceding  measures 
as  indicated  by  equation  (B-V).  These  describing  functions  are 
shown  for  each  of  the  six  stick  configurations  in  Figure  B-10; 
also  shown  in  this  figure  is  the  average  over  four  configurations 
(all  but  heavy  spring  and  heavy  stiff) . Zero  dD  indicates 
1 de<:j/sec  head  rotation  per  g of  platform  acceleration. 

The  curves  shown  in  Figure  D-10  display  appreciably  greater 
frequency  dependency  than  do  the  describing  functions  presented  in 
the  preceding  two  graphs.  The  apparent  cause  of  the  large  swings 
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Figure  B-8.  Deacribina  Functions  Relating  Acceleration  of 
Inboard  Bitebar  to  Platform  Acceleration 
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in  both  amplitude  and  phase  is  that  the  inboard  and  outboard 
accelerometers  are  subjected  to  nearly  the  same  acceleration^ 
for  the  most  part.  As  these  accelerations  approach  identity  in 
amplitude  and  phase,  the  describing  function  will  approach 

zero  amplitude  (i.e.,  a large  negative  dB  value)  and  may  change 
by  180  degrees,  as  can  be  seen  from  inspection  of  equation  (B-7) . 

Equation  (B-8)  was  used  to  compute  the  spectr^-jn  of  the 
acceleration  of  the  head  point-of-regard  with  respect  to  the 
display.  Appropriate  frequency-weighting  and  integration  of 
this  spectrum  provided  estimates  of  relative  rms  displacement 
and  rms  velocity  of  the  head  point-of-regard. 

Estimated  displacement  and  velocity  quantities  for  the  head 
point-of-regard  are  shown  in  Table  B-5a.  As  might  be  expected 
from  the  describing  functions  relating  head  rotation  to  platform 
vibration,  values  predicted  for  the  heavy  spring  and  heavy  stiff 
sticks  were  considerably  greater  than  corresponding  measures 
predicted  for  the  remaining  four  configurations.  Data  provided  by 
the  heavy  spring  and  stiff  sticks  were  considered  anomalous,  as 
explained  below,  and  only  data  from  the  remaining  sticks  was  used 
for  further  analysis. 

The  spectrum  of  the  acceleration  of  the  eye  point-of-regard, 
relative  to  the  display,  was  estimated  by  simply  ignoring  the 
components  of  spectrum  for  head  motion  below  4 Hz.  (That  is,  we 
assumed  the  oculomotor  system  would  compensate  for  head  motion 
below  this  frequency.)  Frequency-weighting  and  integration  were 
then  carried  out  to  obtain  rms  displacement  and  velocity. 

The  results  of  this  analysis,  conducted  for  four  of  the  six  stick 
configurations,  are  shown  in  Table  B-5b.  The  average  rms  dis- 
placement was  normalized  with  respect  to  the  display  gain  to  obtain 
a suitable  visual  noise  pareuneter  for  subsequent  model  analysis. 


Device 

Rms  Motion 

Displacement  i Velocity 
(inches)  { (inches/sec) 

a.  Head  Point* 

-of- Regard 

Light  Spring 

.218 

4.35 

Medium  Spring 

.132 

3.0 

Heavy  Spring 

2.01 

27.4 

Light  Stiff 

.232 

3.73 

Medium  Stiff 

.191 

3.90 

Heavy  Stiff 

1.66 

21.9 

e Point-of-Reqard 


Light  Spring 
Medium  Spring 
Light  Stiff 
Medium  Stiff 
verage 


B . 4 . 3 Investigation  of  Possible  Measurement  Error 


I^s  every  attempt  had  been  made  to  maintain  a constant  bio- 
dynamic environment  throughout  the  experimental  program,  the 
extremely  large  head  motions  (relative  to  the  display)  observed 
for  the  heavy  spring  and  heavy  stiff  stick  configurations  were 
quite  unexpected.  Consequently,  additional  analysis  was  under- 
taken to  determine  whether  or  not  these  results  might  be  the 
rfcs\ilt  of  possible  measurement  error  (where  "measurement  error" 
includes  all  sources  of  error  in  the  experimental,  recording, 
and  analysis  procedures). 

Suspicious  as  to  the  reliability  of  the  data  were  aroused 
by  the  fact  that  data  for  the  heavy  spring  and  heavy  stiff  stick 
configurations  were  obtained  from  the  same  subexperiment  (Experi- 
ment 7),  where  the  stiff  stick  was  at  center  location  and  the 
spring  stick  at  the  side.  Accordingly,  we  proceeded  to  analyze 
the  results  of  the  subsequent  experiment,  in  which  the  same 
sticks  were  used  with  their  locations  reversed.  The  objective 
here  was  to  determine  whether  or  not  there  was  a consistent 
relation  between  head  motion  and  control-stick  design.  If  such 
a consistency  were  to  be  observed  (i.e.,  if  Experiment  8 showed 
the  same  behavior  as  Experiment  7) , we  would  assume  that  the 
subjects  had  changed  their  posture  when  manipulating  these  control 
devices  and  that  such  postural  changes  resulted  in  larger  vibration- 
induced  head  motion.  Alternatively,  if  the  results  of  Experiment  8 
were  similar  to  the  results  of  the  earlier  experiments,  we  would 
question  the  validity  of  the  results  of  Experiment  7. 
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Data  from  Experiment  8 were  not  obtained  in  time  to  ? ' 

digitization  and  analysis  of  the  type  performed  on  preceding 
subexperiments.  Therefore,  a simple  check  was  performed  on  che 
data  of  Experiments  7 and  8 as  well  as  on  data  from  an  earlier 
subexperiment  (Experiment  6) . 

Mean-squared  deviational  scores  were  computed  for  the  inboard 
bitebar  acceleration  and  for  the  difference  between  inboard  and  out- 
board accelerations.  (A  large  difference  score  implies  large  head 
rotations.)  'i.gnals  were  high-pass  filtered  to  eliminate  possible 
d-c  bias  from  the  recorded  signals.  The  primary  variable  of 
interest  was  the  ratio  of  the  difference  score  to  the  score  obtained 
from  the  inboard  accelerometer. 

As  expected  from  the  foregoing  analysis,  the  ratio  score 
obtained  from  the  data  of  Experiment  7 was,  on  the  average,  appre- 
ciably higher  than  the  score  computed  from  data  of  Experiment  6. 

The  data  of  Experiment  8,  however,  confirmed  the  results  of 
Experiment  6 (medium  spring  and  medium  stiff  sticks) . Accordingly, 
we  concluded  that  some  form  of  measurement  error  must  have  been 
present  in  Experiment  7 with  regard  to  the  signals  obtained  from 
either  or  both  of  the  bitebar  accelerometers,  and  these  data  were 
ignored  in  subsequent  analysis. 
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APPENDIX  C 

ANALYSIS  OF  TRACKING  PERFORMANCE 


Tracking  performance  in  both  static  and  vibration  conditions 
is  analyzed  in  some  detail  in  this  Appendix,  and  comparisons  are 
presented  between  experimental  data  and  model  output.  Rms  error 
and  control  scores  are  analy7*?d  first,  followed  by  analysis  of 
frequency-domain  measures. 

Preliminary  analysis  of  performance  scores  was  performed  on 
the  tracking  scores  obtained  from  four  subjects  to  explore  possible 
effects  of  stick  location  on  tracking  performance.  Analysis  of 
variance  revealed  no  statistically  significant  effects.  Accordingly, 
the  bulk  of  the  analysis  of  performance  scores  — and  all  analysis  of 
frequency-domain  measures  — was  performed  on  averages  (where 
possible)  of  measures  obtained  from  the  two  stick  locations. 


C.1  Rms  Performance  Scores 

Rn\s  error  and  control  scores  are  shown  for  each  of  seven 
test  subjects  for  static  and  vibration  tracking  in  Table  C-1. 
Also  shown  are  average  performance  scores  as  well  as  standard 
deviations  and  standard  errors  computed  across  subjects.*  All 
scores  are  in  volts,  referred  to  the  analog  simulator. 


^Definitions  of  statistical  measures  are  given  in  Appendix  A. 


Table  C-1 

Rtns  Performance  Scores  for  Tracking  Variables 


Stick 


Individual  Subjects 


:bh  PH  dm 


Statistics 


DS  [Mean 


a.  Error,  Static  Trackin' 


Light  Spring  .330!  ,3 
Medium  Spring  .33lj  .3 


Heavy  Spring 
Light  Stiff 
Medium  Stiff 
Heavy  Stiff 


.418  .3 

.280  .2 
.262  .244 

.320  .210 


.546  .550  .446  .456  .090  .034 
.483  .355  .497  .410  .083(.031 
.356  .418  .336  .359  .093'. 035 


.319  ,314  .443|.318 
.244  .350  .302  .345 


105 i. 040 
18 


b.  Error,  Vibration  Trackin 


Light  Spring 

Medium  Spring  .435  .522  .679  .502  . 687  I . 652 i . 536 

Heavy  Spring  .508  .424  .548  .403  .651 

Light  Stiff  .471  .495  .647  .480  .545 

Medium  Stiff  .480  .452  .784  .480  .557 

Heavy  Stiff  .438  .417  .513  .426  .485 


c.  Control,  Static  Trackin 


Light  Spring 
Medium  Spring 
Hea\n.  Spring 
Light  Stiff 
Medium  Stiff 
Heavy  Stiff 


.394 

.499 

.390 

.301 

.320 

.391 

.341 

. 397 

.338 

.387 

. 384 
.374 

.377 

.568 

.359 

.408 

.335 
. 345 

.332 
. ^40 

.375 

.376 

.413 

. 350 

.396 

64  .024 
39  .015 
22  .008 
.407  .079  .030 
.382  .026  .010 


d.  Control,  Vibration  Trackin 


Light  Spring 
Medium  Spring 
Heavy  Spring 
Light  Stiff 
Medium  Stiff 
Heavy  Stiff 


.826 

.728  1.071 .4301 .322: .350 
.620  .4761  .529!  .359;  .461 
1.12 
1. 18 

.997!  .7911  .882!  .633 


.528 

r.i54 

1 .568 
.458 
1.C2 

.265 

.102 

.236 

Average  performance  scores  are  shown  in  meaningful  physical 
units  as  well  as  in  volts  in  Table  C-2.  Also  shown  are  the  ratios 
obtained  by  dividing  each  error  and  control  score  obtained  under 
vibration  conditions  by  the  corresponding  performance  measure 
obtained  in  the  static  condition. 

The  ratios  presented  in  Table  C-2c  show  that,  on  balance, 
vibration  caused  proportionately  larger  increases  in  performance 
scores  for  the  stiff-stick  configurations  than  for  the  spring 
sticks.  On  the  average,  spring-stick  error  scores  increased  by 
about  30%,  whereas  stiff-stick  scores  increased  by  about  50%. 

Larger  differences  were  observed  for  the  effects  of  vibration 
on  the  rms  control  score.  Vibration  caused  the  control  score 
to  increase  by  about  30%  on  the  average  for  the  spring  sticks, 
whereas  this  score  increased  by  over  a factor  of  2 for  the  stiff 
sticks.  Particularly  large  increases  were  observed  foi  the  light 
and  medium  stiff  sticks,  where,  as  we  shall  presently  observe, 
vibration  feedthrough  accounted  for  a sizeable  portion  of  the 
control  activity. 

Performance  differences  among  the  six  stick  configurations 
can  best  be  observed  from  Figure  C-1.  Figure  C-la  shows  a somewhat 
lower  rms  error  for  the  stiff  sticks  than  for  the  spring  sticks 
in  the  absence  of  vibration.  Because  of  the  proportionately 
greater  effect  that  vibration  has  for  the  stiff  sticks,  however, 
this  trend  is  essentially  nullified  under  vibration  conditions. 

Figure  C-lb  reveals  a nearly  constant  rms  electrical  control 
signal  for  the  six  sticks  under  static  conditions.  As  noted  above, 
considerable  differences  across  sticks  are  shown  for  the  :^'jntrol 
score  under  vibration  conditions. 
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Table  C-2 


Averaqe  Rms  Performance  Scores 


Variable 


Sprin 


Stick 


a.  Static 


Error  (volts) 
Error  Unches) 
Control  (volts) 
Control  inches) 
Control  (pounds) 


.152 

1.07 


1.74 


Vibration 


Error  (volts) 

Error  Unches) 

Control  (volts) 

Control  Unches) 

Control  (pounds)  i ,422 


1.59 


Error 

Control 


Ratio  of  Vibration  Score  to  Static  Score 


1.39  f 1.25 


Figure  C-lc  shows  that  the  actual  control  force  applied  to 
the  stick  varies  widely  from  stick-to-stick  for  both  static  and 
vibration  tracking.  These  differences,  at  least  for  static 
tracking,  reflect  pilot  adaptation  to  the  electrical  control 
gain.  For  example,  as  the  spring  gradient  is  increased  for  the 
spring  sticks,  the  effective  control  gain  (in  terms  of  volts/pound) 
is  reduced,  and  the  pilot  must  increase  his  rms  control  force  to 
maintain  constant  control  effectiveness. 

Analysis  of  variance  of  the  rms  error  and  (electrical) 
control  scorcw  is  summarized  in  Table  C-3.  Analysis  across  all 
experimental  conditions,  summarized  in  TP.ble  c-3a,  reveals 
svrong  vibration-stick  interactions  for  both  error  and  control 
scores  (significant  at  the  0.001  level).  That  is,  we  can  attach 
statistical  significance  to  the  observation  noted  above;  namely, 
that  the  increase  in  performance  score  due  to  vibration  is  different 
for  the  various  stick  configurations.  Significant  subject-vibration 
and  stick-subject  interactions  are  also  found  from  analysis  of  the 
error  score. 

Because  of  the  significant  interactions,  the  analysis  of 
variance  has  not  for  the  most  part  revealed  strongly  significant 
main  effects.  Differences  in  both  error  and  control  score  due  to 
vibration  are  significant  only  at  the  0.05  level,  and  stick 
configuration  has  a significant  effect  (0.05  level)  only  on  the 
error  score.  These  results  do  not  imply  that  vibration  and  stick 
design  have  no  highly  significant  effects  on  tracking  performance, 
but  that  these  effects  should  be  explored  singly  rather  than  in 
combination. 
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Table  C-3 


Sammary  of  Analysis  of  Variance 


Table  C-3b  sununarizes  the  re&ults  of  analysis  of  variance 
performed  on  the  performance  measures  for  static  tracking  alone 
and  for  vibration  tracking  alone.  We  now  find  that  stick  design 
has  a highly  significant  effect  on  rms  error  in  both  static  and 
vibration  conditions,  whereas  the  effect  on  rms  control  is  signi- 
ficant only  for  vibration  tracking. 

Analysis  of  variance  of  the  data  obtained  from  two  of  the 
stick  configurations  considered  individually  is  summarized  in 
Table  C-3c.  The  effects  of  vibration  on  tracking  error  are  now 
observed  to  be  highly  significant  (.01  level  or  lower)  for  both 
the  spring  and  the  stiff  stick.  A significant  effect  on  rms 
control  score  is  observed  only  for  the  stiff  stick. 

Figure  C-2  shows  the  average  error  and  control  variance 
scores  partitioned  into  input-correlated,  vibration-correlated, 
and  remnant  components,  where  "remnant"  is  defined  as  signal  power 
correlated  with  neither  tracking  nor  vibration  inputs,  vibration- 
correlated  power  contributed  less  than  10%  to  the  total  error 
variance  for  all  stick  configurations.  The  relative  contribution 
to  total  control  variance,  however,  varied  widely  with  the  sticks. 
Vibration-related  power  wa.?  small  for  the  spring  sticks,  but  it 
accounted  for  almost  half  of  the  control  variance  for  the  light 
and  medium  heavy  sticks.  About  25%  of  the  control  variance  was 
vibration-correlated  when  the  heavy  stiff  stick  was  used. 

Differences  in  electrical  control  gain  account  for  the 
different  relative  contributions  of  vibration  feedthrough  observed 
for  the  stiff-stick  configurations.  Because  of  hardware  limitations, 
the  gain  of  the  heavy  stiff  stick  was  about  half  of  that  of  the  other 
two  stiff  sticks.  The  variance  of  the  vibration-correlated  portion 
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of  the  control  power,  measured  as  voltage  output  of  the  control 
device,  was  therefore  reduced  by  & factor  of  nearly  four.  The 
portion  of  the  control  output  related  to  tracking  (including 
remnant)  varied  less  across  the  stiff  sticks,  since  the  pilots 
increased  their  control  forces  to  compensate  for  the  reduced 
control  gain.  Hence,  the  lower  proportional  contribution  of 
feedthrough  seen  for  the  heavy  stiff  stick. 

Differences  in  the  composition  of  control  power  between 
spring  and  stiff  sticks  reflect  differences  in  stick  impedance 
as  well  as  control  gain.  The  reader  is  referred  to  the  main 
text  of  this  report  for  further  discussion  of  the  interaction 
between  stick  parameters  and  performance  measures. 

C.2  Frequency-Domain  Measures 

Frequency-domain  measures  for  each  of  the  six  stick  confi- 
gurations are  shown  in  Figures  C-3  through  C-8  for  tracking  under 
static  and  vibration  conditions.  Each  figure  contains  a pilot 
describing  function  relating  control  output  (in  volts)  to  tracking 
error  (in  volts).  Also  shown  in  each  figure  is  a spectrum  derived 
by  dividing  the  remnant- related  component  of  control  power  (at  a 
given  input  frequency)  by  the  input-correlated  component  of  control 
power.*  This  quantity  provides  an  indication  of  the  pilot's  signal/ 
noise  capabilities  and  is  needed  for  testing  models  for  remnant. 

We  shall  refer  to  this  spectrum  as  the  “ratio  spectrum". 


*the  ratio  is  based  on  what  the  input-correlated  and  remnant  spectra 
would  be  if  the  tracking  input  spectrum  were  continuous  rather  than 
sums-of-sinusoids.  Details  of  this  computation  are  given  in 
Appendix  A. 
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6.  Frequency-Domain  Measures  for  Light  Stiff  Stick 
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7.  Frequency-Domain  Measures  for  Medium  Stiff  Stick 
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. Frequency-Domain  Measures  for  Heavy  Stiff  Stick 
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Model  results  are  compared  witn  experimental  describing 
functions  and  ratio  spectra  in  these  figures.  The  model  results 
do  not  necessarily  represent  the  best  match  to  the  data  obtained 
in  a given  control  situation.  Rather,  we  felt  that  the  results 
of  the  modeling  effort  would  be  more  readily  interpreted  if  an 
acceptable  match  to  the  entire  body  of  data  could  be  obtained 
with  a minimum  set  of  parameter  values.  Therefore,  a single  set 
of  parameter  values  was  found  which  best  matched  the  data  obtained 
from  the  six  stick  configurations  under  static  conditions.  A 
second  set  of  paircimeter  values  was  found  for  the  vibration  data, 
with  the  number  of  parameters  re-evaluated  kept  to  a minimum. 

Three  changes  of  parameter  values  were  needed  to  account  for 
the  effects  of  vibration  on  tracking  performance:  (a)  an  increase 

in  pilot  time  delay  from  0.15  to  0.2  seconds,  (b)  a ten-fold 
increase  in  motor  noise/signal  ratio,  and  (c)  the  introduction  of 
a residual  observation  noise  associated  with  perception  of  tracking 
error.  Values  for  all  pilot-related  model  parameters  are  shown  in 
Table  C-4  for  static  and  vibration  conditions.  Additional  details 
on  the  modeling  effort  — including  the  procedure  for  differentiating 
among  possible  sources  of  vibration-related  interference  — are  given 
in  Appendix  D. 


1 
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Data  from  the  six  control-stick  configurations  show  similar 
changes  due  to  vibration.  The  following  trends  are  observed:  (a) 

the  amplitude  ratio  is  decreased  by  about  2-3  dB  at  the  lowest 
three  input  frequencies  (0.5,  1.25,  and  3 rad/sec),  is  basically 
unaffected  at  6.3  Hz,  and  increases  2-4  dB  at  the  highest  input 
frequency  (10. 5 rad/sec) j phase  lag  is  unaffected  at  lower 
frequencies,  but  shows  an  increase  of  15-40  degrees  at  the 
highest  input  frequency;  the  ratio  spectrum  shows  an  overall 
average  increase  of  about  5 dB  with  no  consistent  change  in 
frequency-dependency.  For  the  most  part,  these  effects  are 
mimicked  quite  well  by  the  model  results. 
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Table  C-4 

Values  for  Pilot-Related  Model  Parameters 


Parameter 


t Function 
e Delay  (s 
or  Noise/S 
ervation  N 
placement 
ocity  Thre 
olacement 


ds) 

} Ratio 
/Signal  Ratio 
shold  (volts) 
d (volts/sec) 
dual  Noise  (vo 
Noise  (volts/ 


Parameter 

Value 

Static 

Vibration 

J = 0^  + .03  ( 
e 

2 — 

3 + qa- 

u •’  u 

0.15 

0.20 

0.004 

0.04 

0.025  (-21dB) 

0.025 (-21dD) 

0.135 

0.135 

0.54 

0.54 

0.0 

0.44 

0.0 

0.0 

*Error  and  control  in  volts#  g chosen  to  provide  a motor  time 
constant  of  0.1  seconds. 


The  describing  function  results  are  summarized  in  terms  of 
gain-crossover  frequency  and  phase  margin  in  Table  C-5.  For  the 
most  part,  gain-crossover  increased  with  increasing  stick  stiffness 
ranging  from  about  3-5  rad/sec  for  static  tracking.  Gain-crossover 
decreased  under  vibration  for  all  stick  configurations,  with  the 
average  decrease  being  about  1.0  rad/sec  for  the  spring  sticks 
and  1.4  rad/sec  for  the  stiff  sticks.  Phase  margin  increased 
under  vibration  by  about  10  degrees  on  the  average. 

T-tests  were  performed  on  the  data  contained  in  Table  C-5 
to  determine  the  statistical  significance  of  vibration-related 
performance  differences.  In  addition,  analysis  of  variance  was 
performed  on  the  ratio  spectrum  at  a single  frequency  (3  rad/sec). 
Table  C-6  shows  that  the  effects  of  vibration  on  gain  crossover, 
phase  margin,  and  the  ratio  of  remnant  to  correlated  control 
power  were  significant  at  the  0.01  level  or  better. 
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^ APPENDIX  D 

^ , MODEL  ANALYSIS 

I 

The  important  results  of  the  model  analysis  have  been 
f summarized  in  preceding  sections  of  this  report.  In  this  Appendix, 

I we  describe  some  of  the  tests  that  were  performed  to  quantify 

I..  certain  model  parameters  and  to  differentiate  among  potential 

I vibration-related  interference  effects. 

r 

I'  D.l  State-Vari abl e Representation  of  the  System 

I 

I A flow  diagram  of  the  control  system  (including  stick-plus- 

|!  hand  dynamics)  is  diagrammed  in  the  main  text  in  Figure  12.  The 

i'  following  identities  are  made  in  deriving  a state-variable 

I,  representation  for  this  system; 


x^  = tracking  input  (a  first-order  noise  process) 

X2  = tracking  error 

x^  = electrical  control  signal 


yj^  = tracking  error 


where  all  state  and  display  variables  are  in  volts. 
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with  the  (theoretical)  tracking  input  assumed  to  be  a 
first-order  noise  process  having  a break  frequency  of  2 rad/sec, 
the  following  set  of  matrix  eauations  defines  system  dynamics: 


-1 

— 

r -■ 

- 

— 

— i 

^1 

-2.0 

0 

0 

0 

^1 

0 

2.0 

2.0 

0 

4.0 

0 

0 

0 

8 

0 

0 

0 

1.0 

* 

^3 

+ 

0 

u+ 

0 

^4 

L.  ~ 

0 

0 

-386.4  K^/M 

-386.4 

B/M 

_ 

386.4 

e'  J 

0 

" 

r 

! 0 

1 

1.0 

0 

0 

<— 1 
X 

V2| 

= 1 0.5 

1 

0 

1.0 

0 

• 

X2 

/3j 

1 

o 

0 

o 

• 

• — 1 

0 

^3 

^4 

w 


where 

w « white  noise  process  with  covariance  equal  to 
the  variance  of  the  input 

K = stick  spring  gradient,  pounds/inch 
s 

B « stick  damping,  pounds/ (inch/sec) 
s 

M = combined  mass  of  stick  and  forearm,  pounds 

The  display  vector  defined  for  this  problem  is  somewhat  non- 
standard in  that  (a)  the  display  quantity  y2  is  equal  to  one-fourth 
the  error  rate,  rather  than  simply  the  error  rate,  and  (b)  the 
pilot's  control  input  is  included  in  the  display  vector.  The 
reason  for  using  one- fourth  the  error  rate  in  this  problem  is 
that  this  quantity  corresponds  directly  to  one  of  the  problem 
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variables  recorued  during  the  analog  simulation.  Model  analysis 
is  unaffected,  except  that  threshold  and  rms  residual  noise  levels 
associated  with  error  rate  should  be  divided  by  4. 

The  third  display  quantity  is  included  to  allow  us  to 
determine  the  pilot  describing  function.  As  presently  implemented, 
the  model  can  predict  describing  functions  between  any  display 
quantity  and  the  "theoretical”  control  variable,  which,  for  this 
problem,  is  the  control  force  "F"  (see  Figure  12.  Therefore, 
we  must  first  predict  the  describing  functions  between  error  and 
force  input  and  between  electrical  input  and  force  input.  We 
then  cascade  these  describing  functions  to  yield  the  desired 
describing  function  relating  electrical  control  input  to  tracking 
error,  which  may  then  bo  compared  with  the  corresponding  measure- 
ment derived  from  the  experimental  data. 

D.2  Static  Tracking 

Model  analysis  was  first  performed  on  the  data  obtained  in 
the  absence  of  vibration.  A single  set  of  model  parameters  was 
derived  to  provide  a good  overall  match  to  the  data  obtained  from 
all  six  control-sti'^'^  configurations.  The  goal  (not  always 
achieved)  was  to  match  the  measured  rms  error  and  control  scores 
to  within  10%  and  to  match  pilot  describing  functions  and  "ratio 
spectra"  (ratio  of  remnant  to  input-rorrelated  stick  power)  with 
1 or  2 dB. 

Determination  of  certain  specific  model  parameters  is 
described  below. 
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D.2.1  Effective  Control  Intertia 

We  first  attempted  to  predict  the  effects  of  stick  design 
on  tracking  performance  using  only  the  stick  mass  as  the  variable 
"M"  in  the  problem  formulation.  This  approach,  however,  did  not 
allow  the  model  to  predict  performance  differences  observed  for 
tho  various  control  configurations.  Much  better  correspondence 
to  experimental  results  was  obtained  when  M was  increased  to 
include  the  m.ass  of  the  forearm  (as  given  in  Allen  et  al.  [7]). 


D.2.2  Cost  Weignting  Coefficients 


Initial  modeling  was  attempted  with  a "cost"  determined  as 
in  previous  analysis  of  single-variable  systems  [1,  2];  namely, 
with  cost  equal  to  a weighted  sum  of  mean-squared  tracking  error 
and  mean-squared  control  rate,  where  the  relative  weighting  was 
chosen  to  yield  a motor  time  constant  of  about  0,1  seconds. 
Reasonably  small  matching  errors  were  found  for  most  of  the  control 
configurations,  but  the  predicted  error  for  the  heavy  spring  stick 
was  about  2 3%  belov/  that  found  experimentally. 


Since  the  control  forces  required  by  this  control  configu- 
ration were  greater  than  those  required  by  the  remaining  five 
configurations,  we  felt  that  the  pilot's  unwillingness  to  generate 
large  control  forces  might  have  been  important  in  this  case. 
Therefore,  a cost  on  control  force  was  added  to  tho  total  cost 
expression,  with  the  associated  weighting  factor  chosen  to  provide 
a good  match  to  the  heavy  spring  data,  (The  weighting  on  control 
rate  was  readjusted  to  maintain  the  motor  time  constant  at  0.1 
seconds.)  Because  of  the  lesser  control  forces  reguired  by  the 
other  control  sticks,  the  match  between  model  output  and  experi- 
mental measurement  was  little  effected  in  these  cases  by  the 
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addition  of  the  cost  weighting  on  control  force.  Accordingly, 
the  augmented  cost  expression  was  used  for  the  remainder  of  the 
model  analysis. 

D.3  Vibration  Tracking 

Having  "calibrated"  pilot-related  model  parameters  for 
static  tracking,  we  then  attempted  to  match  the  data  obtained 
in  vibration  conditions.  Our  goal  was  again  to  obtain  a good 
overall  match  to  all  six  control  sticks  with  a consistent  set  of 
model  parameters  and  to  do  so  with  as  few  parameter  changes  as 
possible. 

D « 3 . 1 Modeling  Visual  Effects 

Having  estimated  relative  rms  displacement  and  velocity  of 
the  eye-point-of- regard  with  regard  to  the  display  (Section  B.4) , 
tests  were  performed  to  determine  the  best  way  of  representing 
the  effects  of  blur  in  terms  of  parameters  of  the  tracking  model. 
These  tests  were  made  for  the  heavy  spring  stick  configuration. 

For  our  first  test,  we  simply  assumed  that  the  effect  of 
blurring  was  to  render  velocity  information  useless.  This 
hypothesis  was  modeled  by  associating  a very  large  noise/signal 
ratio  with  velocity  perception.  Other  model  parameters  remained 
unchanged.  Although  a reasonably  good  match  was  obtained  for 
tracking  error  and  pilot  describing  function,  the  frequency- 
dependency  of  the  ratio  spectrum  was  not  reproduced.  The  model 
predicted  too  little  remnant  power  at  low  frequencies  and  too 
much  at  high  frequencies. 
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For  the  second  test,  we  set  the  rms  residual  noise 
associated  with  displacement  perception  equal  to  the  average 
estimated  rms  eye~display  displacement ; similarly,  residual 
noise  on  rate  perception  was  equated  to  the  estimated  rms 
relative  velocity.  The  results  were  essentially  the  same  as 
in  the  first  test,  with  the  same  gross  mismatch  between  pre- 
dicted and  measured  ratio  spectra. 

Finally,  displacement  residual  noise  was  defined  as  in  the 
second  test,  and  rate  residual  noise  v/as  set  to  zero  (i.e,,  we 
assumed  no  degradation  in  rate  perception) , This  model  change 
alone  produced  an  insufficient  increase  in  tracking  error  and 
control  activity  to  account  for  the  effects  of  vibration.  Good 
results  were  obtained  in  conjunction  with  other  parameter  changes, 
ho'wever,  and  this  treatment  of  visual  effects  was  adopted  for  the 
model  validation  tests  shown  in  the  text  of  this  report. 

D . 3 . 2 Other  Non-Motor  Related  Effects 

Attempts  were  made  to  account  for  the  effects  of  vibration 
through  re-adjustment  of  the  observation  noise/signal  ratio.  This 
model  parameter  is  believed  to  reflect  primarily  central-processing 
limitations,  although  certain  visual  effects  could  also  show  up  as 
a change  in  this  parameter.  Again,  the  model  was  tested  against 
data  from  the  heavy  spring  stick  experiment. 

With  the  remaining  model  parameters  unchanged  from  the 
static  situation,  observation  noise/signal  ratio  was  increased 
to  match  the  rms  error  obtained  during  vibration.  A reasonably 
good  match  was  obtained  to  pilot  gain;  however,  the  predicted 
ratio  spectrum  increased  by  about  2 dB  more  than  was  observed 
experimentally.  We  therefore  concluded  that  changing  observation 
noise  alone  could  not  account  adequately  for  the  effects  of 
vibration  on  tracking  performance. 


i D.3.3  Motor-Related  Effects 

Three  motor-related  model  parameters  were  explored:  (a) 

motor  time  constant,  (b)  time  delay,  and  (c)  motor  noise.  Decreasing 
j the  motor  tire  constant  (i.e.,  increasing  the  weighting  coefficient 
I on  control  rate)  did  not  yield  effects  consistent  with  those 
j observed  experimentally.  For  example,  the  predicted  rms  control 
j score  decreased,  whereas  the  data  showed  about  a 25%  increase  for 
I the  heavy  spring  configuration.  In  addition,  a poor  match  was 
obtained  to  the  frequency-domain  results  at  high  frequencies. 
Consequently,  we  concluded  that  vibration  had  no  significant 
effect  on  this  parameter,  and  the  motor  time  constant  was  hence- 
forth kept  at  0.1  second. 

With  residual  noise  determined  as  described  in  D,3,l  and 
with  other  parameter  values  appropriate  to  static  tracking, 
changes  in  time  delay  and  motor  noise  were  explored.  An 
acceptable  match  to  the  vibration  results  could  not  be  obtained 
through  a change  in  only  one  of  these  parameters.  Although  we 
could  not  precisely  allocate  the  effects  of  vibration  between 
these  two  parameters,  a good  match  to  the  data  was  obtained 
with  a time  delay  of  0.2  seconds  and  a motor  noise/signal  ratio 
of  0.027  (motor  noise  variance  divided  by  the  variance  of  the 
predicted  control  force) . (Corresponding  parameter  values  for 
static  tracking  were  0.15  seconds  and  .0036,  respectively.) 

Model  results  were  then  tested  against  data  obtained  from 
the  remaining  five  control  configurations  to  find  the  set  of 
parameter  values  that  would  provide  the  best  overall  match. 

Residual  noise  and  time  delay  were  quantified  as  described 
above,  and  motor  noise/signal  ratio  was  increased  to  0.04  to 
provide  a better  overall  match.  Remaining  parameter  values  were 
the  same  as  used  in  matching  static  tracking  results. 
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D.3,4  Relation  of  Motor  Noise  to  Control  Variance 


It  was  suggested  in  a recent  analysis  of  lateral-axis 
vibration  that  the  motor  noise  varies  in  fixed  proportion  to  the 
total  control  variance,  where  "total”  variance  is  defined  as  the 
control  variance  related  to  tracking  (including  remnant)  plus 
t''.e  effects  of  vibration  feedthrough  [10]. 

The  spring-stick  data  do  not  appear  to  support  this 
hypothesis.  On  the  average,  vibration  caused  the  total  variance 
to  increase  by  about  a factor  of  2 (i.e. , 3 dB) . Model  analysis 
shows  that  an  increase  in  motor  noise  variance  of  greater  than 
10  dB  is  needed  to  account  for  the  difference  between  vibration 
and  static  tracking.  We  therefore  tentatively  conclude  that, 
while  the  motor  noise  may  vary  with  control  activity,  the  ratio 
of  motor  noise  variance  to  control  variance  is  influenced  by 
the  vibration  environment. 
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APPENDIX  - E ANTHROPOMETRIC  RELATIONSHIPS 

A series  of  15  anthropometric  dimensions  were  measured  on  the  subjects 
used  in  this  study.  These  measurements  were  made  to  determine  the 
comparability  of  the  body-size  distribution  of  the  test  subjects  with 
the  body  size  distribution  of  rated  Air  Force  men  and  to  provide  the 
basis  for  seeking  relationships  between  body-size  characteristics  of  the 
subject  and  effects  of  vibration  on  their  performance.  This  latter 
objective  will  be  pursued  when  a sample  adequate  for  the  required 
analyses  is  available. 

SELECTED  ANTHROPOMETRY  OF  SUBJECTS 
IN  THIS  STUDY  AND  RATED  AIR  FORCE  OFFICERS 


This  Study 

AF  Rated 

Officers  (1967) 

X 

SD 

N 

X 

SD 

N 

Age* 

32.6 

7.2 

8 

29.5 

6.3 

2420 

Weight 

78.6 

7.0 

8 

78.7 

9.7 

2420 

Stature 

175.5 

6.1 

8 

177.3 

6.2 

2420 

Sitting  Height 

91.8 

2.7 

8 

93.2 

3.2 

2420 

Thumb-Tip  Reac’n 

1 80.6 

5.7 

8 

80.3 

4.0 

2420 

*Age  in  years. 

weight  in 

kilograms. 

all 

other  dimensions  in 

centimeti 
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